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Summary— Because of the paucity of data on the phase behaviour of mixtures of compounds 
containing active or polar groups at elevated pressure, the P-V-T-x relations of the iso-butanol- 
n-butanol, methanol n-butanol and diethyl ether n-butanol systems were determined at the 
liquid-vapour phase boundaries from near their atmospheric boiling points to the highest 
temperature and pressure at which the liquid and vapour co-exist. In addition, the vapour 
pressures and saturated liquid and vapour densities of the pure components were determined 
up to their critical points. The data are presented in the form of tables and graphs. 

In so far as is indicated by the data, the phase behaviour of binary systems of the alcohols 
is the same as that exhibited by binary systems of hydrocarbons that belong to a given homologous 
series. 


Résumé—En raison de la pauvreté des données sur les relations entre phases aux pressions 
élevées, pour les mélanges de composés a groupements actifs, ou polaires, les auteurs ont étudié 
les relations -V-T-x pour les phases coexistantes liquide-vapeur depuis le point d’ébullition 
au voisinage de la pression atmosphérique jusqu'aux températures et pressions les plus élevées 
ot les deux phases coexistent, dans les mélanges de butanol normal, avec soit lisobutanol, soit 
le méthanol, enfin loxyde d'éthyle. En outre, on a déterminé jusqu’au point critique, les tensions 
de vapeur et les densité des deux phases pour les quarte constituants purs. Les résultats sont 
exposés sous forme de tableaux et de courbes. 

Dans la limite des faits observés, le comportement des systémes binaires d’alcools reste le 
méme que celui des systémes binaires d’hydrocarbones appartenant a une série homologue 
donnée. 


as exemplified by the OH group in the alcohols 
and the phenols, the COOH group in the acids, 
the NH, and NH groups in the primary and 


The ever increasing use of high pressure in many 
industrial processes has greatly stimulated interest 
in and the need for accurate and extensive 


information on the liquid vapour equilibrium 
relations of mixtures under elevated pressures. 
Up to the present time, the study of these relations 
particularly, in the critical region, has been 
confined, generally, to hydrocarbon systems. 
Such systems, being made of up compounds 
of a distinctly non-polar nature, give no informa- 
tion on the effect of active or polar groups in 
the molecule on these relations. Polar groups, 


secondary amines and other groups, all of which 
contain active-H-atoms, manifest their presence 
in the molecule by the phenomena of association 
and hydration. This phenomenon, under certain 
circumstances due to “hydrogen bonding,” 
materially effects the vapour-liquid composition 
relations in mixtures at atmospheric pressure. 

The present study was undertaken for the pur- 
pose of investigating the vapour-liquid equilibrium 
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relations in mixtures of polar compounds at ele- 
vated pressures. The binary systems, methanol— 
n-butanol, isobutanol-n-but anoland diethyl 
ether—n-butanol were selected in order to study 
respectively, 1) mixtures composed of compounds 
belonging to the same homologous series, 2) 
mixtures of chemical isomers and 3) mixtures of 
physical isomers. The P-T-x relations at the 
liquid-vapour phase boundaries were determined 
from near their atmospheric boiling points to 
highest temperature and pressure at which the 
liquid and vapour phases coexist. 

Hitt and Vax Winkie [5] have determined 
the vapour-liquid equilibrium relations of the 
methanol—n-butanol system at 760 mm pressure. 
So far as is known, none of these systems of 
alcohols have been studied at pressures above 
atmospheric. 


EXPERIMENTAL METHODS 


The methods of determining the vapour pressures, 
the densities and the vapour-liquid equilibrium 
composition relations in the three binary systems, 
iso-butanol—n-butanol, methanol n-butanol and 
diethyl ether—n-butanol were the same as those 
employed and recently described for the propane- 
hydrogen sulphide system [6]. These methods 
involved, very briefly, the determination of the 
P-T curves and temperature — density curves of 
the pure components as well as a series of mixtures 
of known composition in each of the above 
binary systems. The relations between the 
pressure, density, temperature and composition 
at the liquid and vapour phase boundaries 
were then derived by the construction of appro- 
poriate cross plots of these curves. 

The pressure, temperature and density data 
for the construction of the curves were obtained 
by confining a small, air-free sample over mercury 
in a thick-walled glass tube which was surrounded 
by a constant temperature bath. The tube was 
fastened in a mercury-filled compressor block 
with means provided for controlling the pressure 
on the sample. After proper adjustment of the 
pressure and after equilibrium had been attained, 
the volume of the sample was determined by 
measuring the length of the column of sample, 
this length having been related to the volume 


by a prior calibration of the tube. Since the 
phenomena associated with the bubble point 
and dew point are directly observable, the 
pressure, volume and temperature values for the 
saturated states of the sample may be obtained 
as exactly as the observer desires. 


In the present work the temperature was 
measured to within 0.02°C with a copper-con- 
stantan thermocouple and a sensitive potentio- 
meter. For the measurement of the pressure a 
dead weight piston gage similar to that previously 
described [6] was used. Its sensitivity was 0.02 
pound per square inch. The experimental tube 
was constructed of precision bore capillary whose 
volume per millimetre of length amounted to 
0.00320 ce (approximately 2-0mm_ bore). A 
cathetometer, reading to 0-05 mm, was used to 
measure the length of the tube occupied by the 
sample. The methods of calibrating the ther- 
mocouple, pressure gauge and tube have been 
adequately described 


MATERIALS 


High purity commercial samples of methanol, 
n-butanol and isobutanol were used as the 
starting material for the preparation of the pure 
liquid compounds that were studied. Purification 
was effected by distillation in the presence of a 
dessicating agent since the butanols form 
azeotropes with water. For this purpose bright 
magnesium turnings were used with the butanols 
and metallic sodium with methanol. 

The distillation equipment consisted of a 
silvered and vacuum-jacketed column, 25 mm 
in diameter and 2 metres high, filled with glass 
helicies of 3mm, inside diameter and attached 
by a ground joint to a 2-litre boiling flask. The 
still head was provided with a reflux control 
device so that that the reflux ratio could be set and 
accurately measured. All distillations were 
carried out with a reflux ratio of 19 : 1 or higher. 
The receiver in which the distillate was caught 
was attached by a ground joint to the head and 
was equipped with a vent and with a stopcock 
at the bottom for withdrawing the product. 
A drying tube containing phosphorous pentoxide 
dispersed on glass wool was attached to the vent 
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to prevent absorption of moisture from the 
surrounding air. 

Approximately 1500 ml of liquid and 15 grams 
of dessicating agent were charged to the still and 
a middle fraction of about 200ml was taken 
as the purified sample. 

Diethyl ether of “ analytical reagent” grade 
was used without any additional purification 
except that associated with the degassing 
operating which will be described later. 

Prior to their use, a sample of each of the 
purified liquids was tested for comtamination 
by measuring, in the experimental apparatus, 
the difference between the bubble and dew 
point pressures at some elevated temperature. 
The standard boiling point, the density at 0°C 
and the critical temperature, pressure and density 
were also determined and compared with reliable 
published values of these properties. The results 
are summarized in Table I. 


DEGASSING AND PREPARATION 
or MIXTURES 


Air and other non-condensable gases were 
removed from the liquid samples by a series of 
operations which involved freezing with liquid 
nitrogen and pumping off the residual gas over the 
solid, followed by melting and distillation at low 
pressure. Fig. 1 is a diagrammatic sketch of the 
apparatus for the degassing operation as well 
as for the preparation of mixtures of known 
composition in the experimental tube. 

Fifty millilitres of pure liquid were charged to 
bulb 4 and the bulb was attached to the vacuum 


Fig. 1. Apparatus for degassing and preparing mixtures 
of liquid samples. 


line by means of a mercury-sealed ground joint. 
Stopcock 3 was opened, momentarily, in order to 
remove most of the air-rich vapour over the 
liquid, before freezing the sample. As soon 
as the sample was frozen, stopcock 3 was opened 
and the space over the solid evacuated until 
the pressure was less than 1 x 10° mm Hg. 
Cock 3 was then closed and the sample allowed 
to melt, after which 3 was opened, cautiously, 
and vapour pumped off again. This procedure 
was repeated 3 times. Next, the cooling flask 
was brought up around 2, and cocks 5, 7, 9 and 17 
were closed, cocks 1 and 3 opened, and the liquid 
in 4 was distilled and collected in bulb 2, except 
for a small residue which was pumped off into the 
cold trap 18. After evacuating the space over 
the solid, cock 1 was closed, the solid melted and a 
small fraction of the vapour pumped off. Cock 
17 was then closed, 3 opened and the sample 
was distilled and collected in 4, except for a 
small residue which was discarded as before. 
The distillation, back and forth between 4 and 2, 
was repeated three times. The pure degassed 
sample was then stored in bulb 2 in the solid state 
until ready to use. The second liquid component 
was treated in like manner and stored in bulb 6. 

The experimental tube, 16, was attached to 
the vacuum line by a mercury-sealed ground 
joint, 15, as shown in the sketch. 14 was a 
small hair-size capillary which extended into the 
tube to within 8cm of the end. The capillary 
was drawn from the tube of the female half of 
the ground joint, 13. No lubricant was used on 
this joint. Communicating with the capillary 
through side tube, 10, was a mercury reservoir. 
Mercury from the reservoir could be transferred 
to the experimental tube by edjustment of the 
pressure over the mercury in the reservoir. 

In the preparation of a mixture of known 
composition a quantity of liquid was transferred 
by distillation from the storage bulb to the 
small capillary, 8, (0-65 mm bore) whose volume 
per unit length had been previously determined. 
From the length of the tube occupied by the 
sample at 0°C (measured with a cathetometer 
reading to 0-05 mm) and the density of the liquid 
at 0°C, the weight of the sample was calculated. 
A correction was made for the weight of the 
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vapour over the liquid, using the perfect gas 
law and a knowledge of the vapour pressure 
and the total volume of 8 up to cock 7. 

Before transferring the measured sample to 
the experimental tube, the tube was flushed 
several times with a small amount of sample 
which was then discarded to trap 18. By chilling 
the tip of the experimental tube and by manipula- 
tion of the appropriate stopcocks the liquid was 
distilled into 16. In like manner a measured 
quantity of the second liquid was transferred to 16. 
The sample was then allowed to melt and fill 
the end of the tube, after which it was frozen 
again. Mercury was admitted from the fiask 
by admitting air through cock 11. The mercury 
entered the tube through the hair capillary 
just above the frozen sample and rose in the 
tube until the tube was completely filled. In 
this manner any trace of residual gas in the tube 
was displaced by the mercury. 

Cock 9 was closed and the tube and adapter 
12, lowered from the apparatus. Joint 13 was 
disconnected and capillary 14 removed from 
the tube. The tube was then transferred to the 
compressor block and fastened in place, as 
has been described [6]. 


Data 
The P-V-T relations at the saturated states 


n-butanol, diethyl ether—n-butanol, were deter- 
mined. Figs. 2 to 9 are plots of the experimental 
data. The critical temperature and critical 
pressure of the pure compounds and the mixtures 
were observed directly, whereas the pressure, 
temperature and density at the maximum pressure 
point and the maximum temperature point on 
the P-T border curves as well as the critical 
density of the mixtures, were read from large 
scale plots of the data. Critical densities of the 
pure compounds were obtained by the application 
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of the law of rectilinear diameters. The critical 
data for the pure substances and for the mixtures 
are given respectively, in Tables I and II and, 
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Fig. 5. Temperature-density relations in the iso-butanol— 
n-butanol system. 
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Fig. 6. Temperature-density relations in the methanol— 
n-butanol system. 


9 GRITICAL LOCUS 
AN 
= 200 \ 
\ 
mole % 
VA Ne OIETHYL ETHER) 
7 27, 100°00 
29 74-39 
30 47-65 
27-12 
100 
DENSITY g/cm® 


Fig. 7. Temperature-density relations in the ethyl ether— 
n butanol system. 


where possible, are compared with values reperted 
in the literature. Figs. 8 and 9 show graphically 
the relation between the composition and the 
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Fig. 8. Critical pressure-composition relations in binary 
systems of n-butanol with iso-butanol, methanol and 
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Fig. 9. Critical Temperature-composition relations in 
binary systems of n-butanol with iso-butanol, Methanol 
and ethyl ether. 


critical pressure and the critical temperature 
for each of the binary systems. Table III and IV 
(a), (b) and (c) summarize, respectively, the 
saturated properties of the pure compounds 
and the various mixtures at uniform intervals 
of temperature. In Table III selected literature 
values for some of the compounds are included. 

The tables were constructed from the experi- 
mental data by a combination of analytical and 
graphical procedures. The relation between 
pressure and temperature was formulated by 
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Table 1. Physical properties of materials studied. 


i- Butanol Methanol Ethyl ether n- Butanol 
Standard boiling point, °C 
Experimental | 107-83 64-55 34-47 117-73 
Timmermans [9) 107-89 — 64.50 — 34-40 — 117-7 — 
| 108-10 64-75 34-60 118-0 
Density at 0°C, g/cc 
Experimental 08175 0-8100 07363 08245 
TimMEeRMANs (9) 0-8170 08100 — 0-7362 — 08246 — 
08171 0-8101 0-7363 0-8247 
Pressure rise on condensation 
At vapour pressure, Ib./in.* 0-87 0-49 0-94 0-32 
At temperature, °C | 805-7 | 412-6 277-0 246-7 
220-48 153-72 229-70 
Critical pressure, lb./in.? 
Experimental 622-8 1174-4 «529-6 640-4 4 
Literature 1153-6 [9] [10] 523-25 [9] (10) | 711-11 [9) [4] 1955 
| 530-52 [7] 
Critical temperature, °C 
Experimental | 27459 239-43 193-41 289-74 
Literature | 279-07 (9) [1] 240° [10) (9) 193-8 (9) [10] 287-0 [4] 
240-2 [2] 193-9 [7)} [8] 
Critical density, g/cc 
Experimental 0-269 0-272 0-265 0-267 
Literature | 0-275 [9) (2) 
Table 11. Composition and critical constants of mixtures of n-butanol 
| 
Critical point  Marimum temperature point Maximum pressure point 
Composition Temp. P Density | Temp. P i Temp. P Density 
| | 
1786 278-11 6269 
48-97 282-48 632-2 0-268 
23-60 286-25 636-6 0-268 
Mol % MeOH 
83-34 249-84 1133-1 0-275 250-25 1127-6 0-238 249-40 0-315 
69-93 258-22 10775 0-276 25889 1065°6 0-225 257-55 1082-0 0-322 
53-07 267-27 988-1 0-275 268-27 969-4 0-218 266-40 996-5 0-324 
23-81 280-96 798-9 0-271 281-58 785-6 0-222 279-94 806-9 0-315 
Mol % ether 
74-39 219-49 582-8 0-293 220-86 581-8 0-238 220-38 584-0 0-259 
47-65 247-04 625-3 0-284 249-59 622-1 0-224 248-76 625-9 0-254 
27-12 266-04 639-5 0-277 267-41 637-9 0-221 266-78 639-9 0-252 
13-10 278-57 642-4 0-272 279-25 641-1 0-228 279-20 642°5 0-251 
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Fig. 10. Pressure-composition relations. 
Methanol—n-butanol system. 
fitting an equation of the form log P = A — B/T 
to the vapour pressure data for the pure com- 
pounds and to the dew and bubble points for each 
composition investigated ; points in the retrograde 
region were excluded. The equations thus 
obtained (two for each mixture) were used to 
calculated the saturation pressures for each 


g --4 MAXIMUM 
600 PRESSURE 


230 
Z 
4 


20 60 to 
COMPOSITION, mole % ETHYL ETHER —o 


Fig. 11. Pressure-composition relations. 
Ethyl ether—n-butanol. 
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Fig. 12. Y-X diagram for methanol—n-butanol system. 


experimental temperature. The calculated 
pressures were then compared with the measured 
values and graphical residual functions were 
constructed. These residuals could be plotted 
very precisely, since the maximum difference 
between the experimental pressure and that 
calculated from the appropriate equation never 
exceeded 35 pounds per square inch for any 
mixture. The scatter of the residuals about a 
smooth curve seldom exceeded 0-3 pounds 
per square inch and serve as an indication of 
the consistency of the data. 

The densities given in Tables III and IV 
were read from large scale plots of the curves 
shown in Figs. 5, 6 and 7. 
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Table 111. Vapour pressure and orthobaric densities of pure components. 


Diethyl Ether 
Temp. | Pressure Density g/cc 

c | Ib. /in.? | Liquid Vapour 

Exp. Ref. (9) (10) Exp. Ref.(9}(10) Exp. Ref. (9) (10) 

110 118-9 117-61 5996 5942 
120 146-7 145-28 “5822 5764 
130 179-1 177-038 | “5580 
140 216-5 213-69 5435 ‘5885 0449 04488 
150 259-5 256-44 5210 5179 0561 05551 
160 | 308-8 805-10 4954 4947 0702 06911 
170 365-0 361-04 4657 4658 0886 08731 
180 428-9 421-06 4304 4268 1147 1135 
190 501-2 493-34 ‘8687 3663 1690 “1620 

C 193-41 529-6 265 265 
Methanol 
| 
Temp. Pressure Density g/cc 

“Cc lb. sin. Liquid Vapour 

Exp. Ref. (9) (10) Exp. (9) (10) Exp. Ref. (9) [10] 
130 121-5 120-70 | 6718 6770 
140 157-8 156-07 | 6574 6640 
150 201-9 199-86 6420 6495 
160 255-0 251-90 6208 6340 
170 318-3 315-04 6081 6160 
180 393-2 388-46 5891 “5980 0332 03186 
190 481-4 475-98 5681 ‘5770 0412 04010 
200 583-8 575-99 5448 0523 05075 
210 202-7 691-68 “5170 5255 0670 06521 
220 840-8 823-22 4811 0896 08635 
230 1000-4 974-85 4321 4410 ‘1230 “1187 

C 239-43 11744 | 272 272 
n- Butanol 
Temp. Pressure Density g/cc 
| lb. | iqui Vapour 

190 | 112-2 6251 

200 138-9 6109 

210 170-2 5956 

220 | 206-6 ‘5759 

230 248-1 5618 O417 

240 295-4 ‘5422 0507 

250 | 349-0 5197 0625 

260 409-2 4917 0778 

270 477-5 “4585 

2380 | S444 “4146 1817 

C 289-74 | 6W-4 267 267 
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Table III. (cont.) 


iso- Butanol 
Temp. Pressure Density g/cc 
"Cc lb. /in? Liquid Vapour 
220 258-8 5555 0453 
230 308-4 -5338 “0561 
240 364-5 “5032 “0700 
250 427-9 “4789 “0880 
260 500°1 “4412 11438 
270 582-2 2830 “1614 
C 27459 622-8 | 269 “269 


C = critical point. 


Table 1V. Summary of pressure, temperature and density data at phase boundaries. 
(a) iso-Butanol—n-Butanol Systems 


Temp. Saturated Liquid | Saturated Vapour 


Composition 
% i-Butanol | °C P density | P density 
| | Ib. /in? g/ec Ib. /in.? g/cc 
77-86 | 200 | 169-0 5963 
210 205-2 -5795 | 
220 246-8 244-2 0418 
230 294-1 54038 291-6 0513 
240 347-8 845-4 0638 
| 250 408-5 406-3 0796 
| 260 427-5 . 475-4 1022 
270 555-6 4096 553-5 -1389 
| 278-11 626-9 269 626-9 
48-97 | 190 128-6 6158 
| 200 158-4 . 
210 192-8 
220 232-2 
230 | 277-2 5485 273-9 0472 
240 328-3 -5261 325-2 0585 
250 | 286-4 323-6 0731 
260 452-2 47038 449-9 0924 
270 526-1 524-6 
280 609-5 -3600 608-8 1825 
C 282-48 632-2 268 632-2 268 
23-60 200 148-1 6054 
210 180-9 5898 
220 218-7 -5781 
230 262-0 -5550 259-8 0438 
240 310-9 “5847 308-7 0541 
250 366-5 364-4 0671 
260 429-7 -4803 427-5 0841 
270 500-8 4442 498-6 -1079 
280 580-9 -3910 578-6 1434 
C 286-25 636-6 268 636-6 268 


C = critical point. 
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Table 1V-b (cont.) 
(b) Methanol—n-Butanol System 


Composition Saturated Liquid Saturated Vapour 
mol % MeOH P density P density 
lb. gyec lb. g/cc 
53-07 150 126-3 6713 
160 159-5 6589 
170 198-8 6458 
180 245-0 6312 
190 299-0 6157 
200 361-2 +5983 
210 432-4 -5789 
220 513-4 377-5 0451 
230 604-7 “5300 456-6 0570 
240 704-7 4992 | 550-3 0720 
250 814-5 -4600 662-5 0923 
260 932-0 4048 796-2 “1258 
MP 266-4 996.5 324 VOL. 
C 267-27 988-1 275 988-1 275 4 
MT 268-3 969-4 218 1955 
23-81 170 126-4 6514 
180 157-3 6378 
190 193-3 6232 
200 235-0 GOTT 
210 282-8 “5903 
220 337-3 “5715 
230 398-8 5504 312-4 0451 
240 466-5 -5273 374-4 0563 
250 540-7 446-0 
260 623-1 4695 528-6 0897 
270 713-9 623-8 “1188 
MP 279-9 806-9 B15 
280-96 798-6 27 798-6 271 
MT 281-6 785-6 222 
83-34 140 136-6 6735 
150 1743 
160 219-9 
170 274-3 
180 338-5 6138 
190 413-6 -5957 
200 500:5 “5751 416-9 
210 600-6 5517 511-2 0510 
220 715-2 -5230 622-6 0652 
230 844-9 +4882 753-7 0855 
240 989-0 908-6 “1169 
MP 240-4 1185-1 B15 
C 249-84 1183-1 275 1183-1 275 
MT 250-3 1127-6 “238 
MP = Maximum Pressure Point. 
MT = Maximum Temperature Point. 
C = Critical Point. 
10 
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(b) Methanol—n-Butanol System 


Composition Temp. Saturated Liquid | Saturated Vapour 
mol % MeOH °c P density P density 
lb. /in.? g/cc lb. /in.2 
69-93 150 152-0 6664 
160 191-9 “6525 
170 239-3 6372 
180 295-0 4210 
190 360-2 -6030 
200 435-7 5831 
210 522-2 5611 398°8 0413 
220 620-2 -5367 486-9 “0527 
230 730-1 591-3 
240 852-7 -4698 714-4 -0876 
250 985-6 “4154 864-9 1213 
MP 257-6 1082-0 322 
C 258-22 1077-5 1077-51 
MT 258-9 1065-6 225 


Table 1V-c (cont.) 
(c) Ethyl Ether —n-Butanol System 


74:39 130 141-7 6569 
140 171-1 6393 
150 204-9 -6209 
160 243-5 -6006 
170 287-2 209-5 0419 
180 336-6 5528 258-1 0529 
190 391°8 “5214 815°5 0673 
200 453-1 383-1 “C878 
210 520-9 4340 463-0 1194 
C 219-49 582-8 582-8 
MP 220-4 584-0 259 
MT 220-9 581-8 238 
47-65 140 125-0 
150 149°8 6460 
160 178-5 
170 211-2 “6122 
180 247-6 
190 289-9 “5727 
200 336-7 5491 237-2 -0419 
210 389-2 -5220 289-7 0531 
220 447-1 -4908 351-3 0677 
230 510-6 -4530 423:2 -0886 
240 579-9 *B992 507-4 1201 
C 247-94 625°3 284 625-3 284 
MP 248-8 625-9 254 
MT 249-6 622-1 224 


MP = Maximum Pressure Point. 
MT = Maximum Temperature Point. 
C = Critical Point. 
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Table 1V-c (cont.) 
(c) Ethyl Ether—n-Butanol System 


Composition Temp. Saturated Liquid | Saturated Vapour 
mol % Ether °c P density P density 
Ib. /in®. g/ee lb. /in?. g/vv 
27-12 160 131-1 6487 | 
170 | 157-0 6327 
180 186-7 6162 
190 220-7 “5983 
200 259-1 
210 202-6 227-2 0398 
220 351-3 5354 274-8 -0480 
230 | 405-5 5101 329-9 0602 
| 240 465-4 -4802 393-3 0762 
250 531-4 “4431 466-4 
260 603-6 ‘8873 552-6 1840 
C 266-04 629-5 ‘277 639-5 ‘277 
MP 266-8 639-9 252 
MT 267-4 637-9 221 
13-10 180 128-5 6259 
| 190 166-6 6102 
200 199-1 “5942 
210 236-2 5771 
220 278-4 “5587 233-8 
230 326-1 “5379 280-9 0471 
240 379-8 “5132 335-0 0588 
250 439-8 -4834 396-7 0746 
260 506-5 “4489 467-1 0953 
270 580-3 4014 547-0 +1273 
C 278-58 642-4 ‘272 642-4 272 
| MP 279-2 642-5 251 
| MT 279-3 641-1 


MP = Maximum Pressure Point. 
MT = Maximum Temperature Point. 
C = Critical Point. 


Figs. 10 and 11 are constant temperature plots of 
the relation between pressure and composition at 
the bubble and dew points of the methanol— 
n-butanol and diethyl ether—n-butanol systems. 

Table V (a), (b) and (c) presents the liquid- 
vapour composition data for the three systems, 
in terms of the familiar y, z and K (K = y/z) 
parameters. Figs. 12, 13, 14 and 15 show some 
of these data graphically. The table was con- 
structed from values read directly from large 
seale plots of temperature-composition isobars. 

In general, only the last significant figure of 


any value appearing in the tables is subject to 
uncertainty. The pressure and temperature 
values are estimated to be accurate to within 
+ 0-5 pounds per square inch and + 0-05°C, 
respectively. The uncertainty in the density 
values is somewhere between 0-0020 gram per 
ce for the liquid and 0.00020 gram per cc for the 
vapour, though some of the recorded values in 
the critical region of the mixtures may be subject 
to an error as large as 0.006 gram per cc. The 
compositions given for the mixtures are accurate 
to within 0-05 mole per cent. 


VOL. 
4 
1955 


Liquid-vapour equilibria in the iso-butanol—n-butanol, methanol— n-butanol and diethyl ether — n-butanol systems 


Table V. Temperature-composition relation at constant Table V 
pressure. (c) Ethyl Ether—n-Butanol System 
(a) iso-Butanol—n- Butanol 
P | Temp z y K einer \Kn-puon 
b.jin® 
P Temp. z 
Djin §=—°C 200 | 1358 | 1000 | 1000 | 1-00 
140 90-6 97-0 1-07 319 
400 245-75 100-0 100-0 1-00 150 71-9 90-0 1-25 -356 
247 90-2 91-7 | 1-02 160 56-4 818 | 1-45 417 
240 746 | 170 43-2 724 | 1-68 486 
251 58-7 62-6 1-07 906 180 31-5 61-5 1-95 
253 43:1 46 | 100 085 190 21-6 489 | 2-26 652 
255 27-4 02 | 110 “061 200 12-9 343 | 2-66 134 
257 18 182 | 112 006 210 5-3 168 | 3-17 879 
258.5 00 218-3 0-0 0-0 1-00 
600 27203 | 100-0 100-0 1-00 300 158-3 100-0 100-0 1-00 
274 85-2 864 | 101 919 160 96-2 98-7 1-03 342 
276 70:3 72:3 1-03 933 170 778 91-0 1-17 405 
L. 278 55-4 579 | 1-05 os 180 63-0 822 | 1-30 481 
: 280 40-5 42-9 1-06 960 190 49-8 723 | 1-45 552 
55 282 25-7 273 | (106 978 200 37-7 61-3 | 1-63 621 
284 10-7 113 | (1-06 993 210 26-8 48-9 1-82 698 
285-4 0-0 0-0 1-00 220 1g} 347 2-03 788 
230 8-3 189 | 228 884 
240-9 0-0 0-0 1-00 
1800 
£ #00 | 175-7 100-0 100-0 1-00 
g q 180 92-4 96-1 1-04 518 
100¢ 190 76-8 86-7 1-13 566 
200 62-3 77-0 1-24 610 
\ \ 210 49-2 66-6 1-35 657 
600 \ \ Qk 220 37-2 553 | 1-49 ‘712 
q \ 270 230 26-2 42-8 1-63 ‘775 
240 16-2 28-7 1-77 ‘851 
NX Nese \ 250 70 137 | 1-96 928 
\ 258-6 0-0 0-0 1-00 
230 ‘\230 °C 
q \ 500 | 189-8 1000 | 100-0 1-00 
200 200 84-9 99-0 1-06 “662 
210 210 70-7 79:8 1-13 689 
\ i 220 57-4 692 | 121 726 
ne 490 230 45-0 58-0 | 1-29 764 
100 240 33-3 46-1 1-38 ‘808 
Kmeon = x 260 12-0 198 | 1-65 911 
EQUILIBRIUM RATIO —e 270 2-6 +8 1-85 977 
Fig. 14. Pressure-equilibrium ratio curves at constant 378-1 
temperature. Methanol—n-butanol system. 600 230-0 65-1 65-1 1-00 1-00 
240 30-9 56-8 1-12 880 
Discussion 250 38-2 45-4 1-19 883 
The combination of the vapour pressure curves 260 266 28-5 ed pen 
of the pure components and the P-T border — bo + — 
280 5-3 7-3 1-42 977 
curves of the mixtures form a space diagram 285-4 0-0 0-0 1-00 


representing the P-T-x relations of the system. 
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Table ¥ 
(b) Methanol—n-Butanol System 
P Temp | zy | Kyou 
lb. /in. 
200 | 1496 | 100-0 100-0 1-00 
160 736 89-1 1-21 412 
170 | 53-6 77-9 1-45 476 
180 88-0 65-7 1-73 553 
190 25-7 522 | 203 643 
200 15-7 36-1 2-30 758 
| 210 6-9 172 | 249 ‘871 
| 215 2-7 69 256 956 
218-3 0-0 0-0 100 
400 | 180-8 100-0 100-0 1-00 
| 190 80-0 91-2 1440 
620 B12 1-31 495 
47-5 70-0 1-47 571 
220 35-2 57-7 1-64 653 
230 242 44-0 1-82 739 
| 240 148 29-2 1-97 “831 
250 6-6 13-8 2-09 923 
| 955 27 | 58 | 215 968 
258-6 oo 00 1-00 
600 | 201-4 100-0 100-0 100 
210 886 91-8 1-10 500 
220 66-7 81-6 1-22 553 
52-6 70-4 1-34 624 
40-9 58-8 14 697 
250 80-5 46-5 | 1-52 ‘770 
211 338 839 
280 42 1-71 970 
285-41 oo 1-00 
| | 
soo | «2172 1000 100-0 | 1-00 
2300 780 870125 
687 763 1-20 653 
| 0 518 65-3 1-27 713 
2006 «40-9 537 | 131 -788 
280 | 22-9 25-1 | 1-10 971 
282-3 214 21-4 1-00 1-00 
1000 | 230-0 | 100-0 100-0 100 
240 84-3 89-1 1-06 694 
250 10 «77-9 1-10 762 
260 508 (65-7 1-10 853 
267-7 | 540 540 1-00 | 1-00 


These diagrams for the diethyl ether-—n-butanol, 
methanol — n-butanol and iso-butanol_- n-butanol 


systems are shown in Figs. 2, 3 and 4 respectively. 
The broken line in each of the diagrams is the 
critical locus or envelope curve defining the critical 
state for all mixtures of the components. The 
shape of the critical locus appears to be a function, 
principally, of the relative volatility of the two 
pure components so long as they do not form an 
azeotrope. If it be assumed that the critical 
temperature of any mixture is the mole fraction 
average of the critical temperatures of the com- 
ponents, then, on the assumption that the critical 
pressure of the mixture is the sum of the partial 
pressures of the components at the critical 
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15. Pressure-equilibrium ratio curves at constant 
temperature. Ethyl ether—-n-butanol system. 


Fig. 


temperature, the shape of the critical locus curve 
can be predicted, qualitatively. Thus, if the 
relative volatility of the components is close to 
unity, that is, if the vapour pressure curves lie 
close together and are parallel, and their critical 
pressures are about the same, the critical locus 
would be predicted to be a straight line. These 
conditions are fulfilled and the prediction amply 
verified by the results on systems composed of 
chemical isomers such as the iso-butanol— 
n-butanol system, Fig. 4. If the vapour pressure 
curves of the components are not too widely 
separated and the component having the lower 
critical temperature has the higher critical 
pressure, then the critical locus will be convex 
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upwards as in the case of the methanol— 
n-butanol system, Fig. 3. If the relative volatility 
of the components were greater, the locus would 
be expected to contain a point of maximum 
pressure. A similar result would be expected if 
the component having the greater critical tem- 
perature has the greater critical pressure. In 
both cases there must be a minimum value of the 
relative volatility at which the maximum pressure 
point just appears. The diethyl ether—n-butanol 
system illustrates such a case, since the point of 
maximum pressure is only 2-5 pounds per square 
inch higher than the critical pressure of n-butanol. 

The critical loci of the three binary systems 
were calculated in the manner indicated, the 
“vapor pressures”’ of iso-butanol, methanol 
and diethyl ether, above their critical tempera- 
tures being estimated by extrapolating the vapor 
pressure curve plotted as the logarithm of the 
pressure versus the reciprocal of the absolute 
temperature. The calculated curves were found 
to be in surprisingly good agreement with the 
experimental curves considering the approximate 
nature of the assumptions. The maximum 
deviation in the critical pressure amounted to 
approximately 5% for the diethyl ether—n- 
butanol and methanol—n-butanol systems and 
to about 1% in the case of the iso-butanol— 
n-butanol system. This method of calculation, 
however, is not one of general application. In 
hydrocarbon systems, for example, the deviations 
between calculated and experimental critical 
pressures may be very large. 

Information on the magnitude of the deviation 
of the three systems from ideal behaviour is 
given by a comparison of the vapour-liquid 
equilibrium ratio, K = y/z, for n-butanol in 
each of the systems, with the value calculated 
[3] from ideal solution laws using the fugacities 
instead of vapour pressures. The data are 
presented in Fig. 16 as curves of log P versus 
log K. As will be noted, at the lower pressures 
all of the curves lie very close together indicating 
an approach to ideal behaviour. At higher 
pressures, however, the nature of the associated 
component exerts a major influence and the 
experimental curves pass through a minimum 
K value and terminate at a pressure where K = 1, 
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EQUILIBRIUM RATIO ——> 
Fig. 16. Comparison of equilibrium ratio curves 
for n-butanol in different liquids: i—iso-butanol, 
M—methanol, E—ethy] ether, F—ideal solution. 


which is the critical pressure for a specific mixture 
of the two components. At the critical point, 
all mixtures deviate from ideal solutions in their 
behaviour. In the systems concerned, n-butanol 
shows the least deviation in mixtures with 
iso-butanol and the greatest deviation in mixtures 
with ethyl ether. 

The difference between the values of the vapour 
pressure and saturated densities for methanol 
and diethyl ether, and those selected from the 
literature by TimmMerRMANs [10] as the most 
reliable (Table III), is worthy of comment. 
The data selected are those of Youne [11]. 
For the determination of the pressure, Younc 
used an air manometer which was calibrated at a 
relatively low pressure [12]. The densities were 
determined in the same manner as in the present 
work, except that the length of the calibrated 
tube which was occupied by the sample, was 
measured on a millimetre scale etched on the tube. 
The same method was employed to obtain the 
volume of the air in the air manometer. These 
methods of measuring pressure and length are 
not as precise as those employed in this work 
and the difference in precision is believed to 
account for the difference in values of the vapour 
pressure and densities recorded in Table III. 
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Summary—The amount of SO, absorbed by individual drops of water and of several hydro- 
carbons has been measured continuously during their formation at the tip of a capillary tube. 
The times of formation varied between 1 and 40 seconds. A simplified calculation using the 
* penetration "’ theory explains the experimetal results fairly well. At long formation times, 
however, the rate of absorption is somewhat increased by free convection currents, whereas at 
short times of formation a considerable increase in absorption rate may be found as a consequence 
of the mixing produced by the jet of liquid leaving the capillary. A tentative criterion for the 
occurrence of the latter effect is given. 


Résumé—La quantité de SO, absorbée par les gouttes individuelles d'eau et de plusieurs hydro- 
carbures a été mesurée d'une facon continue pendant leur formation a l’extrémité d'un tube 
eapillaire. La période de formation de la goutte a été variée entre 1 et 40 secondes. Un calcul 
simplifié basé sur la théorie de la “ pénétration * explique assez bien les résultats expérimentaux. 
Cependent, pour de longues périodes de formation, la vitesse d’absorption est augmentée iégére- 
ment par des courants libres de convection, tandis que 4 des périodes courtes de formation on 
peut trouver une augmentation considérable de la vitesse d’absorption par suite du mélange 
produit par le jet de liquide quittant le capillaire. On propose un critére déterminant la naissance 


de ce dernier effet. 


INTRODUCTION 

The rate of mass transfer between a drop and a 
surrounding fluid has been subject to numerous 
experimental studies. In most of these experi- 
ments an extraction or absorption process has 
been carried out in a small-scale spray column. 
In such a process a distinction can be made 
between the exchange of material during the 
formation of the drop, during its free movement 
through the continuous phase and during coales- 
cence. Generally, however, only the overall result 
of the exchange is measured. The sum of the 
beginning and end effects can then be estimated 
by extrapolation to zero drop life from results 
for the various times during which the free drop 
is in contact with its surroundings. 

This procedure has been followed by Suer- 
woop [1], Licur and Conway [2], and West et 
al. [8] for liquid-liquid extraction. Their results 
show that the rate of extraction during drop 
formation is relatively rapid. Quantitatively, 


their data do not agree very well and they are 
contradictory as to the influence of the rate of 
formation of the drops on the initial extraction. 

CouLson and SKINNER [4] have used a different 
method for measuring the initial extraction of a 
drop. They formed their drops on the tip of a 
jet, suspending them in the continuous phase, and 
let them disappear back into the same jet 
immediately after formation. They found smaller 
initial extraction rates than the previous workers. 
The same conclusion was reached by Licut and 
Pansinc [5]. They showed that the extrapolation 
of extraction measurements according to SHER- 
woop’s [1] method may give values which are 
higher than the actual amount extracted during 
the growth of the drop; the release and subse- 
quent vibration of the drop may account for an 
appreciable portion of the entire initial extraction 
effect. 

Absorption of CO, by water drops was investi- 
gated by Wuirman [6] in a small column of 
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constant height. He varied the time of formation 
and by extrapolation to zero formation time 
found the amount of gas absorbed during the 
free fall of the drop. Assuming this to be constant 
he calculated the absorption during formation by 
subtracting this constant value from the total 
amount absorbed. Drxon and Russexu [7], on 
the other hand, found that during the free fall 
of drops the rate of absorption of gas increased 
when the rate of formation of the drops was 
increased. They measured the initial effect by 
making the height of fall as small as possible, and 
found that under conditions of strong turbulence 
in the growing drop quite high absorption rates 
could occur. 

The present investigation was undertaken to 
obtain a better understanding of the internal 
resistance of a drop to mass transfer during its 
formation. “An experimental method was 
developed by which the amount of gas absorbed 
could be measured continuously during this 
period. To this end, drops were formed at the 
tip of a capillary in an enclosed volume containing 
pure SO,. The amount of gas absorbed was 
determined by the measurement of pressure or 
volume changes. The results were compared with 
simple diffusion theory, which, for the sake of 
clarity, will be treated in the next section. 


ABSORPTION BY A Drop 


It is well known that the rate of diffusion of a 
solute into a semi-infinite, stagnant medium which 
has a concentration ¢c = c, at t = 0 and a con- 
centration c* at its boundary for t > 0, is given 


by . 
N” = (c* — 


Here N” represents the rate of diffusion through 
the boundary per unit surface and D the coefficient 
of diffusion of the solute. This equation is also 
valid for diffusion into a layer of finite thickness, 
as long as the depth of penetration of the solute 
is relatively small. From this starting point, 
which has previously been used by Hicsre [8] 
and Danckwerts [9] for treating gas absorption 
phenomena, it is possible to make a rough estimate 
of the absorption of gas by a drop being formed 


(1) 


at the tip of a capillary if the surrounding gas is 
not subject to diffusional resistance. 

The surface area A of the drop increases with 
time. Let the area increment in the time interval 
between @ and @ + d@ be dA (@), then at the 
time ¢(> @) the rate of mass transfer through 
this area is : 


dN = (c* J d4(0), (2) 


and the total rate of absorption at time ¢ : 


AW 
dA (0) 


N (t) = (c* — /— 
4, Ve — 6) 


where A (o) = 0. 

In this derivation it has been assumed that the 
depth of penetration is small compared with the 
radius of curvature of the surface, that there is 
no turbulence within the drop, and that the new 
surface being formed is really “ fresh” surface 
with concentration c,. The latter assumption is, 
of course, questionable, and entirely disregards 
possible diffusion parallel to the liquid surface. 

For the case in which the drop is being formed 
at constant flow rate we have, for the volume of 
the drop : 


(3) 


V= (4) 
The relationship between V and the surface area 
A of the drop depends essentially on the condi- 
tions and the stage of formation. Drxon and 
Russe [7] measured the volumes and areas of 
drops on capillary tips (2-8 mm O.D.), and found 
a reasonably constant ratio of volume to surface 
area. Since in the present investigation similar 
capillaries were used, it was considered justified 
to put 
A=V/r, (5) 

r, being a dimensional constant which was found 
to be 0-61 mm from Dixon and Russexv’s data. 

After substitution from (4) and (5), equation 
(3) can be integrated to give : 


2 
N (t) = (et — (6) 
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The amount of gas, m (t), absorbed after a time t 
then follows from : 


 & 


Consequently, the degree of saturation after a 
time ¢ is : 


EXPERIMENTAL Part 


As can be seen from Fig. 1, the drops were 
formed at the end of a thick-walled glass capillary 
at a constant rate by displacing the liquid with 
mercury. In its turn the mercury was displaced 
by a plunger, which was pushed into the mercury 
by a mechanical drive. The flow rates could be 
adjusted by varying the thickness of the plunger 
and by using different plunger velocities. The 
volume of the drop at any stage of its formation 
could be derived from readings of the mercury 
level in the calibrated vertical tube (diam. about 


Fig. 1. 


This relationship cannot be expected to be 
valid for relative saturation values greater than 
about 0-2, because the penetration of the solute 
then becomes too great and the assumption of a 
semi-infinite medium is no longer valid. In 
addition to this, the other assumptions have been 
such that equations (7) and (8) can be expected 
to have only qualitative significance. 


Apparatus used for formation times longer than 5 seconds. 


1 mm). SO, was used throughout the experiments 
as the gas to be absorbed. Only two capillaries 
were used : 


Tip No. 1, inner diameter 0-85 mm, 
outer diameter 2-6 mm 

Tip No. 2, inner diameter 0-35 mm, 
outer diameter 2-6 mm 
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Fig. 2. Apparatus used for formation times shorter than 5 seconds. 


The water drops fell off when their mass had 
become about 40 mg whereas the maximum drop 
mass of the hydrocarbons used was about 16 mg. 
Care was taken that the liquid did not wet the 
outer wall of the capillary ; this was inevitable, 
however, with heptane and in a few cases with 
nonane, thereby causing much higher rates of 
absorption. Such results have not been reported 
here. 

For low rates of formation (40 mg in more than 
5 seconds) the vessel in which the absorption took 
place was connected with a manometer filled with 
acetone (Fig. 1). At the beginning of an experi- 
ment the absorption vessel, whose volume was 
accurately known, was entirely filled with SO, 
gas which was saturated with respect to the 
liquid used. By measuring the pressure change 
occurring during the growth of the drop, the 
amount of SO, absorbed could be calculated after 
correction for the drop volume. The pressure 
changes were always less than 10 cm of acetone. 

For these slow runs the response time of the 
manometer proved to be negligible. When the 
drops were formed in less than 5 seconds, the 
amount of gas absorbed was determined by 
measuring the change of volume at constant 
pressure, which was indicated by the movement 
of a drop of acetone in a calibrated glass capillary 


(Fig. 2). As the liquid experiments took 1 to 5 
seconds photographs were taken of the position 
of the acetone drop and the level of the mercury 
column at intervals of 1/16th second. A syn- 
chronous clock making one revolution in 5 seconds 
was included in the picture. 

All measurements were carried out in triplicate, 
and only average results are given here. In the 
slow runs the reproducibility of m (t) proved to 
be well within 10%. For the fastest runs the 
spread in the results was larger than 10% in a 
few cases. The volumes absorbed were between 
0-4 and 0.005 cm® of gas. 

During a run the gas and liquid temperatures 
were kept constant (in general at 20°C), but as a 
result of the absorption the drops were heated 
about 2 degrees centigrade above the ambient 
temperature. In view of the uncertainties involved 
in the theory and in the diffusion coefficients for 
the liquids, no attempt was made to take this 
temperature rise into consideration. 


Resutts Water Drops 


Fig. 3 shows a typical set of results for the weight 
of SO, absorbed, m (t), as a function of the time 
during which the drop was growing. It appears 
that the rate of absorption is higher for the 
shorter times of formation, but that the total 


| i © 
ai 
ia 
1955 


——> t pec] 


Fig. 3. Typical curves of the amount absorbed vs. time 
during drop formation ; water. 


amount absorbed is less than for the longer times 
of formation. 

These data have been plotted in Fig. 4 as 
m (t) vs. Vt}, in accordance with equation (7). The 
ratio of m(t) to Vt! is very nearly constant, so 
that from these measurements a proportionality 
factor K can be derived : 


K. (9) 


According to the simple theory which has been 
derived above, K should be 


(9a) 


and should thus be independent of the time of 
formation. From Fig. 4 it is seen, however, that 
K decreases for shorter times of formation. 

The values of K,, obtained for absorption by 
water with zero initial SO, concentration, have 
been assembled in Table 1. This table shows that 
with slow formation there is no significant differ- 
ence between the values of K, obtained with the 
wide and with the narrow capillary. 

In both cases a decrease in K is observed as 
the formation time becomes less. At high flow 
rates, however, the drop formed on tip No. 2 
absorbs much faster than the drop on tip No. 1. 
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Fig. 4. Data of Fig. 3 plotted vs. Vi according to equation 
(7). 


Table 1. Absorption of SO, in growing water drops at 20°C. 


= 0. 
Ky, 10° g/cm sect Approx. total 
formation time 
10-8 cm3/sec | Tip No. 1| Tip No.2 sec 
10 51 5-4 40 
15 51 5-3 27 
2-0 48 45 20 
2-9 41 46 14 
40 3-9 41 10 
8-0 3-8 3-7 5 
3-4 42 24 
33-6 3:5 51 1-2 


This effect is caused by the greater circulation in 
the drop on tip 2, due to the jet of liquid leaving 
the narrow capillary. 

The increase of K, with decreasing values of 
®, can be ascribed to internal circulation in the 
drop in the reverse direction, which is due to 
density gradients. The density of a saturated 
solution of SO, is 7% greater than that of water 
which is sufficient to cause appreciable convection 
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currents in these small drops. This was demon- 
strated by separate experiments in which phenol- 
pthalein was added to water which had been 
made slightly alkaline. Since the water was 
decolorized upon absorbing the SO, the flow 
pattern within the drop was made visible. 

Fig. 5a shows an intermediate stage of absorp- 
tion by a water drop hanging in SO, gas (time of 
formation about 20 seconds). It is evident that 
there is sufficient time during the formation for 
the liquid containing SO, to flow to the bottom 
of the drop. A similar photograph is shown in 
Fig. 5b; here a (rising) drop was formed in CCl, 
which contained SO,. It is rather surprising that 
with this convection taking place the propor- 
tionality between m (t) and Vt}, derived from the 
very simple theoretical model given before, still 
obtains. 

The intensity of free convection could be 
reduced by decreasing the difference in density 
between the surface and the bulk of the liquid 
by using a partially saturated solution of SO, in 
water. The results are indicated in Table 2. From 
the last column it is seen that about 10°, of the 
absorption rate at long formation times and 
Cy = 0 can be attributed to free convection. 


Table 2. Absorption of SO, in growing water drops having 
different initial SO, concentrations at 20°C; 4, = 1:5 x 
10-3 cm® /sec., Tip No. 1. 


} 
K | 3 
Cy/C* 10° g /em® sect | 10-3 g /em® sect 
| 
0 5-0 (Ky) 30 
0-17 40 48 
0-38 29 £7 
2-5 | 45 
0-72 13 | +6 
10 0 


In order to check equation (9a) the following 
values were used: For c* at 20°C the value of 
0.0995 g/cm® was taken as indicated by Jonn- 
sTONE and Leppia [10] for the equilibrium 
solubility of SO, in water. Thus it was assumed 
that the hydration and subsequent partial 
ionization of SO, in water are relatively fast. For 


D, the literature indicates values between 1-8 and 
1-5 x 10° cm?/sec. Measurements by the authors 
(see also footnote to Table 4) on diffusion of SO, 
into a stagnant water surface gave as a result : 
c* 1/D = 0.87 x 10% g/cm? sec! + 5%, yielding 
D = 1-4 x 10% cm? /sec + 10%. With these 
values, r, = 0-061cm and C, = 0, K, is found 
to be according to equation (9a) : 


K, = 46 x 10° g/cm! sec}, 


Thus the theory, in spite of its rough-and-ready 
character, provides a reasonable agreement with 
the K,-values of Table 1. 

From the dependence of c* on temperature and 
the assumtiopn that D is proportional to T'/y, it 
can be calculated that K, should be increased by 
nearly 2 per cent on lowering the temperature by 
1° in the region between 10° and 20°C. This was 
confirmed by a number of measurements at 10°, 
15° and 20°C, the results of which are given in 
Table 3. 

A few experiments in which a small amount of 
surface active agent was added to the water 
showed that this had no significant influence on 
the value of K. In these runs the final drop 
volume was of course smaller, owing to the lower 
surface tension. 


ABSORPTION BY HyDROCARBON Drops 


Normal alkanes were used in order to vary the 
physical properties of the absorption system. 
Moreover, in these substances SO, is only 
physically dissolved. Some properties of these 
liquids are given in Table 4. 

The measurements of m (t) as a function of drop 
life gave results similar to those for water. In 
most cases m (t) proved to be exactly proportional 
to Vt. Only at high rates of formation with tip 
No. 2 did a plot of the data yield curves slightly 
convex to the Vit-axis. This can be ascribed to 
the fact that the jet of liquid from the capillary 
promotes the absorption more effectively in the 
initial stages than in the final stages of the drop 
formation. The least viscous liquid, n-heptane, 
behaved very exceptionally, giving a much faster 
absorption than was expected (up to 10 times 
the calculated rate for the short formation times). 
Unlike the other liquids, its surface was seen to 
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Fig. 5. Occurrence of free convection in absorption of SO, 

in slowly formed water drops : (a) drop hanging in SO, 

atmosphere, formation time about 20 sec.; (b) drop rising 
in CCl, containing SO,, formation time about 10 sec. 
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Table 3. Absorption of SOz in growing water drops at 10°, 15° and 20°C, co = 0, Tip No. 1. 


Ky (10°) | Kg(15°) | Ko (20°) K, (10°)/Kg (20°) Kg (15°)/Kg (20°) 
10-8 cm /sec x 108 x 108 x 108 meas. _ calc. meas. _ calc. 

10 5-52 5-16 1-07 

1-5 5-76 5-73 4-98 1-16 1-15 

2-0 ~ 4-75 461 103 f 

29 4-74 4°34 4-00 1-19 1-09 


Table 4. Properties of hydrocarbons (20°C). 


Solubility Diffusivity 
Substance Viscosity of SO,* of SO,** 

g/emsec| 10-5 cm? /sec 
n-heptane 0-41 0-0288 2-7 
n-nonane 0-71 0-0265 2-5 
n-decane 0-92 0-0250 24 
n-dodecane 1-51 0-0226 20 
n-hexadecane 3-34 0-0188 14 


* These values have been taken from earlier work done 
by the first author in the Koninklijke /Shell-Laboratorium, 
Amsterdam. 

** These values were determined experimentally from 
the rate of absorption of SO, exposed to a flat, fresh inter- 
face of the stagnant liquid. The equipment for measuring 
the rate of absorption was similar to that shown in Fig. 2. 


The diffusivities are thought to have an accuracy of = 


+ 10%. 
Table 5. Absorption of SO, in growing hydrocarbon drops at 
20°C, tg = 0. 
Ko, 10° g/cm? sect 
10°% Tip 
em®/sec | No. | nonane | decane | dodecane 
10 1 20 16 1-2 0-92 
1-5 1 20 1-6 1-4 0-98 
2-0 1 19 14 11 0-89 
29 1 18 14 ll 0-90 
40 1 15 14 11 0-89 
8-0 1 1-4 1-4 1-2 0-94 
16-8 1 13 1-1 10 0-70 
16-8 2 29 13 10 0-70 
33-6 1 18 14 0-85 0-85 
33-6 2 6-3 2-7 14 0-85 
K, calculated, 
eq. (9a) 164 | 154 | 124 0-88 


be in a state of high agitation during absorption. 
This effect was independent of the rate of forma- 
tion. 

Table 5 and Fig. 6, which summarize the 
results for the hydrocarbons, do not include the 
measurements with heptane. 

The values of K, for the hydrocarbons show the 
same trends as those observed for water: an 
increase with smaller values of ®, (due to free 
convection) and with higher values of ®, (due to 
jet mixing, especially noticeable with tip No. 2) 
with a minimum at intermediate rates of forma- 
tion. Both effects tend to diminish as the viscosity 
increases. For hexadecane K, is practically con- 
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Fig. 6. Review of results in terms of Ky as a function of 
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stant. There is a rough agreement between the 
calculated and the measured values of Ky. 


GENERAL Discussion 


All experimental results can be correlated accord- 
ing to equation (8) by plotting the relative satura- 
tion (¢ — ¢,)/(e* — ¢,) versus (Dt/r2)*. This has 
been done in Fig. 7, where the points indicate 
the degree of saturation for the moment at which 


the observed phenomena. The effect of mixing 
by free convection and jet action can be seen to 
increase the relative saturation at long and at 
short formation times respectively, for the least 
viscous liquids. According to Dixon and 
RussE.L’s measurements, ih the case of water r, 
does not vary much with the outer diameter of 
the tip in the range between 0-2 and 0-8 cm. Thus 
equation (8) is believed to give a reasonable 
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Fig. 7. Correlation of all results according to equation (8). 


the drop has a volume of 0-020 cm*. Thus each 
point is representative of an experimental run 
which would give a straight line passing through 
that point and the origin. The choice of 0-020 cm* 
is arbitrary ; it is about half the final volume of 
the water drops as they are released and about 
equal to the final volume of the hydrocarbon 
drops. 

It appears that equation (8) roughly describes 


estimate of the amount absorbed by a drop 
during its formation when the liquid side resist- 
ance is predominant and as long as the jet mixing 
effect is unimportant. The following criterion is 
proposed for the point at which the jet effect 
becomes important : the force exerted by the jet 
on the liquid in the drop is counteracted by viscous 
shear and by free convection due to density 


differences occurring during absorption. 
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If one neglects this effect of convection, one 
comes to the conclusion that the average velocity 
of the liquid leaving the capillary, a, the viscosity 
of the liquid, yu, the inner diameter of the capillary, 
d,, and the drop size must be the most important 
factors affecting the jet action. In these measure- 
ments u varied from 0-2 to 85 cm/sec and p» from 
0-7 to 88cP; d; was either 0-085 or 0-035 cm. 
From the reported data it was found that for a 
value of 


~ 49 to 50 

the rate of absorption begins to increase as a 
result of jet action. Whether there is also an 
influence of the average drop size cannot be 
determined from these data since the outer 
diameter of the capillary tube was not changed. 

Dixon and [7] expressed the “degree 
of turbulence” in their growing drops in terms 
of 1/d; at constant formation time (roughly 
constant ®,). This is not in contradiction to the 
dimensionless criterion mentioned above. By 
analysing their reported data one can find that 
above a value of pu d;/yz of about 50 the rate 
of CO, absorption increases when the formation 
time is kept constant. However, the relative 
saturation found by them for “quietly”? formed 
drops is about twice the value calculated by 
means of equation (8). This is probably due to 


the fact that in their results the absorption during 
release of the drop is included, whereas in the 
measurements reported here only absorption 
during growth was measured. 
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NOTATION 


A = surface area (cm?) 
¢ = concentration of SO, in liquid (g/cm*) 
¢ = mean concentration 
Co = initial concentration 
c* = saturation concentration 
d,; = inner diameter of capillary tube 
D = coefficient of diffusion of SO, in liquid (em*/sec) 
K = constant defined by equation (9); Ky for cg = 0 
m(t) = amount of SO, absorbed after ¢ sec (g) 
N” = diffusion rate (g/cm? sec) 
N = absorption rate (g/sec) 
1, = V/A for growing drops, here 0-061 cm 
t = time (sec) 
T = absolute temperature (°K) 
u = average linear velocity (cm/sec) 
V = volume of liquid drop (em*) 
= dynamic viscosity (g/cm sec) 
p = liquid density (g/em*) 
@ = variable time (sec) 
®, = volumetric rate of liquid flow (em*/sec) 
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The use of a sedimentation cell in the sampling of aerosols 
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Summary—The limited range of conditions under which a reliable sample is obtained are 
discussed, and recommendations are made for the design of the cell. 


Résumé—Les auteurs discutent |'étendue limitée de conditions dans lesquelles on peut obtenir 
un échantillon valable et donnent des conseils pour le projet de la cellule. 


1. INTRODUCTION 


The sedimentation cell has long been considered 
a reliable means of sampling an aerosol and has 
been used as a standard in the appraisal of other 
sampling devices [1], [2]. It is generally con- 
structed so that a volume of aerosol can be 
enclosed and then allowed to settle on to a glass 
slide at the bottom of the cell. The particles can 
later be examined and counted with the aid of a 
microscope. 

The cell used in the present work is shown 
diagrammatically in Fig. 1 and was similar to 


that described by WuytLaw-Gray and Patrer- 
son [3]. It consists of a spacer P with a circular 
orifice C, and two slides S, and S, which are moved 
simultaneously to isolate a sample of aerosol; a 
glass plate is fitted in the lower slide S,. 


The cell was found to suffer from certain 


* Now with the Safety in Mines Research Establishment. 


limitations which had not previously been 
reported. 

2. RANGE OF APPLICATION 
The various sampling errors which can occur 
during use of the instrument will now be discussed 
and the range of conditions over which it gives a 
representative sample will be indicated. 

2.1. When particles settle at random on a 
horizontal surface, a certain fraction « overlap 
one another. Armirace [4] has shown that when 
mn, identical circular discs of radius a fall at 
random on unit area : 


= 1 — a® + termsin(n, a*)* ete. (1) 


Thus if the particles are deposited from a column 
of height d of aerosol of numerical concentration n 


— nd (2) 


provided only a small proportion of the area is 
covered by particles. 

If «, is the largest permissible fraction of 
particles overlapping, the depth of the cell must 
be not greater than : 


d = (1 —«,)/27a*n (3) 


2.2. The sample on the slide will not be 
representative if an appreciable proportion of the 
aerosol settles as the cell is closed. 

In Fig. 2, a particle of radius a and settling 
velocity v is shown initially near the leading edge 
of the upper surface of the cell. If the slide is 
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closed with a velocity u, the particle strikes the 
lower surface at a distance x from its leading 


edge, where : 


r=ud/v (4) 
- 
— | 
Fig. 2 


and the distribution of particles over the slide 
is not uniform. If the diameter of the cell is 
2r, the time of closing is 2r/u and an element 
of the slide at a distance x from the leading edge 
is exposed for a time (27 — z)/u. The number of 
particles deposited per unit area at 2 is then 
nvu(2r —z)/u, and at the instant of closure 
n(d —av/u) particles per unit area are still to 
be deposited. Hence the total deposition N per 
unit area at x is: 


N=n[d +2(r —2)v/u] (5) 


for values of less than r + ud /2r. At greater 
distances there is no deposition. 

If the aerosol is polydisperse, n and v are 
functions of the radius a of the particles, and the 
number deposited per unit area at x is : 


N = n(a) [d + 2(r —2)v(a)/u] (6) 


where n(a) is the numerical concentration of 
particles of radius a and v (a) is their sedimenta- 
tion velocity. An example of the distribution 
obtained in a sample of an oil mist is shown in 
Fig. 3. 

The numerical concentration n in the aerosol 
can be determined graphically from the inter- 
cepts on the axis in Fig. 3 since : 


n = N,(1 —2r/2,)/d (7) 


where N, is the intercept on the N axis and 2, is 
that on the 2 axis. 


27 


2.3. Although it would appear from equation 
(4) that the rate of closing u should be as high 
as possible, excessively rapid movement of the 
slides can cause turbulence with the result that 
particles are deposited on the walls and the upper 
slide, giving an uneven distribution as shown in 
Fig. 3. 


PARTICLES OF RADII < hip 


RELATIVE NUMBERS OF PARTICLES 


° 100 200 
DISTANCE FROM LEADING EDGE OF SLIDE 


Fig. 3 


2.4. The sample will not be representative if 
an appreciable proportion of the aerosol diffuses 
to the walls of the cell. As the diameter of the 
cell is arranged to be very much greater than its 
depth, only the top and bottom surfaces need be 
considered. With material of density 1 gram per 
cu. cm in a shallow cell, this loss becomes signi- 
ficant with particles less than about 0-5 microns 
in radius. 

The depth d of the cell should in any case be 
so chosen that it is at least ten times as great 
as the root mean square displacement of the 
particle (y*)# due to its Brownian motion during 


the sedimentation period, t. Since =f and 


y? = 2 Dt, from Ernsrein’s equation [5], the 
minimum depth for the cell is given by 


d = 200 — = 200 — (8) 
v mg 
ince D = — 
67an 67a7n 
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where & is Boltzmann’s constant 
T is absolute temperature 
g is the acceleration due to gravity 
m is the mass of the particle 
a is the radius of the particle 
7 is the viscosity of the gas. 


2.5. When the aerosol is composed of liquid 
droplets, they may spread. If a droplet is treated 
as a plano-convex lens and its focal length is 
measured, the true radius may be calculated by 
a method due to May [6], provided the radius is 
more than about 5 microns. 


8. DeEsIGN oF THE CELL 


An estimate must first be made of the concen- 
tration of the aerosol and equation (3) used to 
calculate the maximum depth for which over- 
lapping is inappreciable. If this is greater than 


the minimum depth given by equation (8) the 
method of sampling is feasible, i.e., if : 


(1 — €,)/27 > 200 kT /mg (9) 


To prevent loss of larger particles near the leading 
edge, the diameter of the cell should be greater 
than ud/v (from equation 4). For the purpose 
of calculation, the speed of closing must be 
estimated. A large diameter is advantageous 
since the loss of particles to the side walls is 
reduced. 
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On the steady state fractionation of multicomponent and complex 
mixtures in an ideal cascade 


Part 1—Analytic solution of the equations for general mixtures 
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Summary—The equations for the steady state rectification of mixtures which obey a rather 
general vapour-liquid equilibrium law in an ideal cascade are solved in closed form. Various 
kinds of mixtures are examined. From the mathematical standpoint the problem presents some 
interesting features, first, because it may be linearized and, second, it can be reduced to an 
eigenvalue problem which in turn can be solved by a new transform, the transform being finite or 
infinite depending upon the number of components in the mixture. Methods for the numerical 
evaluation of the formulae as well as applications to practice will be reported elsewhere. The 
calculation of the minimum reflux will be presented in Part II. 


Résumé—Dans le cas d'une cascade idéale, auteur a résolu les équations de la rectification 
en état stationnaire de mélanges qui obéissent @ une loi d’équilibre liquide-vapeur suffisamment 
générale. Il examine diverses sortes de mélanges. Du point de vue mathématique, le probléme 
présente quelques caractéristiques intéressantes : (I°) parce qu'il peut étre mis sous forme 
linéaire et (2°) parce qu'il peut étre ramené & un probléme “ de valeurs particuliéres * qui peut 
étre résolu, & son tour, par une nouvelle transformation, cette transformation étant finie ou 
infinie, suivant le nombre de composants du mélange. L’auteur donne-ra ailleurs des méthodes 
pour l’évaluation numérique des formules aussi bien que des applications pratiques. I présentera 
le calcul du reflux minimum dans la deuxiéme partie. 


employed to solve the basic equations for the ideal 
column are both rather unusual and sometimes 
even intriguing. Methods for the numerical 
evaluation of the formulae will be reported 
elsewhere. 


1. INTRODUCTION 


To be able to design a fractionating column, or 
cascade, and to predict its behaviour at steady 
state has been a well known problem which is 
still far from solved. Considerable progress to- 
wards the eventual solution of this problem has 
been made during the past years and many of 
the factors which influence the operation of 
cascades are beginning to become clearer and to 
be understood. It still remains true, nevertheless, 
that the physical model itself is too involved to 


2. Tue CoLuMnN. 
Tue Basic EQUATIONS FOR THE 
RECTIFICATION OF DiscrRETE MIXTURES 


A simple cascade can be represented schematically 
by Fig. 1. Vapour travels up the column and 


be treated exactly and it is for this reason that 
several mathematical models have been con- 
structed and used with some success. It is the 
purpose of this paper to exhibit solutions of the 
most commonly used model, the ideal column. 
From a mathematical standpoint the problem 
has some interesting features because the methods 


liquid moves in the downward direction the two 
being contacted in stages so that the two streams 
leaving a stage are in equilibrium, more or less. 
The column can be divided into two sections, the 
upper part or enriching section and the lower part, 
or stripping section, the two sections being 
separated by the introduction of the feed. For 
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Fig. 1. Schematic diagram of a distillation cascade. 


the moment, only discrete mixtures, that is, 
mixtures with a finite number of distinguishable 
components, will be considered. By a mass 
balance around the enevlope E, the following 
equation is obtained, 


(2.1) 


= ayli) + 


Vass 
where 


Yn (i) = mole fraction of the ith component in 
the vapour stream leaving the (n+ 1)th 
stage. 

z,(i) = mole fraction of the ith component in 
the liquid stream leaving the nth 
stage. 

L,, = moles of liquid leaving nth stage. 

V.., = moles of vapour leaving (n + 1)th stage. 

D = moles of overhead stream. 

r(i) = mole fraction of ith component in the 

overhead stream. 


It will be assumed that the streams leaving a 
single stage are in thermodynamic equilibrium 


and that the equilibrium is described by the 
equation* 


Palt) Yn?) 
where p, (i) is related to the vapour pressure of 
the component, the vapour-liquid equilibrium 
constant, or the relative volatility, and m denotes 
the total number of components in the mixture. 
Combining Eqs. (2.1) and (2.2) there results 


z, (i) = (2.2) 


where 
Ve, =L, + D 
L, (2.4) 


It will be seen from the above equations that 
there are two independent variables, namely, i 
and n, where i is the component designating 
variable and n denotes the stage number. Cer- 
tainly i and n are integers such that 1 s i s m 
and n = 0, n = 0 denoting the condensing stage. 

Now, in general, Eq. (2.3) is misleading to the 
extent that R, and p,(i) are functions not only 
of n but also of the temperature and composition 
of the liquid mixture in each stage. To take this 
into consideration in an analytical way presents 
insurmountable difficulties so that the model will 
be simplified further: viz., it will be assumed 
that p, (i) = p(i) for all n and R, = R for all 
n. The basic equation, then, which must be 
solved is given by 

R Phi) yal) 


and it is desired to obtain an analytic expression 

for y,(i) subject to the condition that 
y,(i) = 

a self evident relation, ~ ‘ich is valid whenever 

a total condenser is used. This model, although 

somewhat restrictive in nature, holds for a wide 


class of mixtures such as isotopes and hydro- 
carbons. Its use in many cases is of doubtful 


+ (2.5) 


(2.6) 


* Summations are over the components in the mixture 


from 1 to m. JZ, denotes that an obvious term must be 
omitted from the summation. 
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value for precise results, but, as in the case of all 
analytic formulae, the qualitative and asymptotic 
deductions are important. A detailed discussion 
of cascades and distillation processes may be 
found in a wide variety of references. The solution 
of Eq. (2.5) will now be found for discrete mix- 
tures. The solution for non-discrete mixtures will 
be found later in the paper. 


8. Tue SoLuTION OF THE PROBLEM FOR 
Discrete MixTuREsS 

1. The general solution of the equations for the 

enriching section. Our first task is to linearize 

Kq. (2.5). If one considers the equation for n=1, 

then it follows that 

R pli) 


R+1Zp(ijafi) R+1 

If in this formula alone p(i) is replaced by Y,(i) 
and y,(i) is substituted into the expression for 
y,(?) there results 


yt) = 


P(t) = 
[R pli) ¥ (i) + pli) (i) 

P(A) ysl?) = 
agli) [R pli) + pli) ¥ 
R+1 


1 agi) ¥,(i) 


~ R+12 
if Y,(i) is defined to be 
Y,(i) = R pli) ¥,(i) + pli) Y,fi) 2,(i) 
Hence it is clear that if one defines 
= pli) RY, (i) + pli) £¥,(i) agli) 


then 
Pi) = EY, ali) ali) 


It is convenient for later work to have Y,(i) = 
RX,(i) so that 


1X, (i) aol) 


P(t) = $1 Lagi) X,.(i) (3.1) 


and 
= pli) + pli) (8-2) 
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It is to be noted here that Eq. (3.2) is an 
equation of partial differences to which must be 
appended the condition 


xo = (3.8) 


Thus the solution of the non-linear difference 
equation, Eq. (2.5), may be solved once the solu- 
tion to the linear difference equation, Eq. (3.2) has 
been found. Equation (8.2) may be solved since 
it may be written 


X,.,(i) = H X,(i) (3.4) 


where H is a linear operator. The iterative method 
applies and the solution is 


X,(i) = pli) (3.5) 


an expression which is rather cumbersome to 
apply. This iterative solution is identical with 
what is commonly called the “ plate by plate ” 
method. 

If one assumes the variables in Eq. (3.2) may 
be separated then X,(i) = v(i) N(n) and 


v(i) N (n + 1) = R pli) ofi) N(n) + 


P(i) N(n) xo(i) 
and 
N(n + 1) _ Rpli)v(i) + pli) Z 
v(i) 


By the usual argument 


N(n + 1) = RA N(n) 
and 


H = R phi) oli) + pi) E = 


where RA has been chosen as the representation 
for the eigenvalue. Therefore if it is assumed, 
for the moment, that there are g eigenvalues then 
the solution may be written in the form 


X,(i) -£ C, (RA,)** (3.6) 
<1 


v,(i) being the representation for the Ath eigen- 
function. Also since 


— 
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it follows that if this formula is multiplied by 
z,(i) and summed on i from 1 to m, then 


PAi) Toft) 
R. (3.7) 
This equation determines the eigenvalues and 
since p(i) is assumed to be a single valued function 
of i for the discrete case, there are m, g = m, 
distinct eigenvalues all of which are real and 
positive as shown by the intersections on Fig. 2. 


A-py) 


Fig. 2. Qualitative plot of equation for determination of 
eigen values. 


It is clear that v,(i) may be defined by 


as) 
— Pli) (3.8) 
From (3.7) and (8.8) it follows that 
(v,, = 5," (3.9) 


where 5,*’ is the Kronecker symbol and the inner 
product (v,, v,-) is defined by 


(vy, = v,(j) 

Thus the eigenvalues are non-degenerate and 

the eigenfunctions form an orthogonal set. That 

the set is complete can be readily proved, for if 

Sli) is an arbitrary function of i defined at 

i = 1, 2,--m then it can be expanded in terms 
of the eigenfunctions v,(i), i.e., 


(3.10) 


v, 
kel 


for 1 sism. The above is but a set of m 


linearly independent linear equations in m 
unknowns and so possesses a unique solution in 
the ¢,’s since the eigenfunctions are linearly 
independent. In order to determine the C,’s in 
Eq. (3.6) it need be noted only that 


Hence 


C, = (v,, (3.11) 


It can easily be shown from Eqs. (3.1) and (3.2) 
that 


2, (i) = ali) X.i) (3.12) 
and therefore 
C, 0, (i) 
z,(i) = ke (3.13) 
Cy 
k=l 
Eq. (3.13) represents the solution of the problem 
originally posed. 


Thus far the one to one correspondence of the 
symbol i and the component in the mixture has 
been left arbitrary. For convenience the p(i) will 
be ordered so that 


p(1) > p(2) > > > p(m). 


The eigenvalues will also be ordered in such a 
way that 
p(k) < A, < p(k — 1) 


and from Fig. 2 this is clearly possible. 

2. The stripping section. The equations for 
the stripping section are essentially identical with 
those of the enriching section. It can readily be 
shown that 


R*—1Z R*—1 
where L* = R*W and 2,*(i) is the mole fraction 
of the ith component in the liquid mixture leaving 
the zeroth plate (reboiler) at the bottom of the 
cascade, For the sake of generality it has been 
assumed that p() is different in the two sections 
of the column, p(i) + p*(i). Eq. (8.14) can be 
solved in the same manner as (2.5) and it is 
found that 
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> C,° oy s- 


To*(i) 


m 
C,* 
k=l 


1 v,*(i) 


z,*(i) = (3.15) 
where 

pri) 
P*(t) — de 
and the ¢,’s are the m roots of 


v,*(i) (3.16) 


_ 
p*(j) — 
ordered in such a way that 


p*(k + 1) < d& < p*(k) 


(3.17) 


and 
C,° 
j= 1 (PX) — 
Finally at the plate where the raw feed is 


introduced, where n = N and s = S, one must 
have 


(3.18) 


ry(i) = 
or by combining (3.13) and (3.15) 


Tolt) 
R 
kel 


kel 


If x,(i) denotes the mole fraction of the ith 
component of the feed, then by a mass balance 
around the whole cascade (per mole of feed) 


xi) = Dali) + W 2*,(i) (3.20) 


From Eqs. (3.19) and (3.20) it is apparent that 
one cannot specify a priori either r,(i) or 2*,(i). As 
a matter of fact, if R, the composition of the feed 
and its thermal condition are given and if it is 
stated, moreover, that the feed is to be introduced 
at such a plate in the column as to make the total 
number of plates for a given separation a mini- 
mum, then it is easy to show that there are two 
degrees of freedom remaining, irrespective of the 
number of components in the mixture. These are 


normally exhausted by specifying z,(i) at i = 8 
and z*,(i) at i = 

3. A useful asymptotic expression. A useful 
asymptotic form for z,({i) can be obtained if 
is made.small. First the function 


Pj) tl 


is defined. Now if w(i) > 0, which is generally 
the only case of interest, and if z{7) is small, then 


(3.21) 


A; = pli) (1 (3.22) 
provided that 
j w(i) 
or, more precisely, that 
Also 
ej) _ Pj) tJ) 
[Ai — PU)? ~ pli) rol?) 
__ 


and if (3.23) is satisfied, then from the above and 
(3.11) it follows that 


= P(t) oft) 8.24 

(wt) 

provided that 2,(i) is sufficiently small. In addi- 
tion, if 

(3.25) 


aoli) Cy vi) < < wli) C, AN" 
kei 
then one easily obtains from the above and (3.13) 


Equation (3.26) is an asymptotic expression which 
is valid only if z9(i) is sufficiently small so that 
Eqs. (3.23) and (8.25) are satisfied, and if w(i)>0. 
Naturally, an entirely analogous asymptotic for- 
mula holds for x,*(i) if z9*(i) is sufficiently small. 
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Thus far the separation of mixtures with a 
finite number of components in an ideal cascade 
has been examined and the basic equation (2.5) 
has been solved analytically without recourse to 
approximations. Similar problems have been 
discussed by Harpert [1], Unperwoop [2], 
Murpocn [8], Hirayama [4], and the present 
authors [5], by different methods. It is felt that 
the present method is more straightforward and 
has the particular advantage that it may be 
extended to more complex mixtures and in 
particular to mixtures whicli have a continuous 
distribution function. 


4. Tue SOLUTION OF THE PROBLEM FOR 
CompLex MIXTURES 


1. Continuous Miztures. In recent years the 
theoretical treatment of the separation of mix- 
tures with an indefinite number of components, 
that is, where the concentration of the various 
components can be approximated by a distribu- 
tion curve, has received some attention. Bow- 
MAN [6] was the first to transform the equations 
which describe a state for a discrete number of 
components so that they could equally well apply 
for the continuous case; he then used those 
equations to study some systems, namely, batch 
distillation and flash vapourization [7]. It is the 
purpose of the latter part of this paper to present 
Bowman's equations for the rectification of 
mixtures with an indefinite number of components 
and to solve those equations both by the iterative 
method and by a scheme analogous to that 
presented earlier and involving the use of a new 
integral transform. The theoretical treatment, 
aside from its mathematical interest, will be found 
useful in studying the rectification of petroleum 
mixtures and, in many instances, accurate enough 
to be used in the discrete case where the number 
of components is so large that the equations 
developed earlier would be somewhat cumbersome 
to use. 

It will be supposed that each component in the 
continuous mixture is present in infinitesimal 
amount. For discrete mixtures the component 
designating variable was the integer i in p(i). For 
continuous mixtures it is natural to employ the 
quantity p itself as the variable designating the 


component since this is, in general, unique at a 
fixed temperature and pressure. It is apparent 
that the mole fraction cannot be used to specify 
the composition of a mixture but it is much more 
convenient to specify the composition of the 
mixture by a distribution function, defined as 
follows : Let z,(p) dp be the mole fraction in the 
liquid phase leaving the nth stage of material 
having p values between p and p+ dp. The 
vapour distribution function is defined ana- 
logously. Methods for calculating the distribution 
function of a given mixture from its experimentally 
obtained true boiling curve have been exhaust- 
ively discussed by Harsert [8] and Bowman [9]. 


The conditions 
(4.1) 
oO 
(4.2) 


are self evident. Further conditions to be imposed 
on 2,(p) and y,(p) are: (a) single valued, (b) 
integrable, and (c) lim [p*2,(p)]}=0 and 


pro 
lim [p* y,(p)] = 0 for all e. 
pro 
Equations for complex mixtures analogous to 


those for discrete mixtures for the enriching 
section are 


P ap 
and 


R+1 
Jp dp 


where R and 2,(p), the distribution function of 
the overhead product, have been specified a 
priori. The linearization of (4.4) follows in the 
same manner as previously. A function X,(p) is 
defined such that 

X,(P) (4.5) 


P ~Ri12 
dp 
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where 


X,4:(P) = RpX,(p) + p J rp) X,(p)dp (4.6) 


or symbolically 

X,+:(p) = X,(p) (4.7) 
where H is a linear operator and 

X\(p) = p/ R. (4.8) 
The iteration of (4.7) gives 

X,(p) = (4.9) 


This form of the solution although not very useful 
here does settle the question concerning existence 
and uniqueness of the solution, for if (4.9) is 
written out in detail 


= dp + Clo) 


dp + -- + | rol p) dp 


where the C,(p) are polynomials in p and the 
integrals are the successive moments of the over- 
head distribution function. Certainly for machine 
calculations these integrals are not difficult to 
compute. These integrals will exist from our 
previous hypothesis ec. 

However, an integral transform can be defined 
which will generate a solution which has some 
advantages. If f(i) is a function of i which is 
defined for integral i, 1 s i S m, then a finite 
transform T [f(i)] = f(k) may be defined by 

— — tel?) 

T [f(i)] = fk) — pli) (4.10) 
where 29(i) is specified and A, is one of the roots 
of (3.7). Since the v,(i) form a complete orthogonal 
set 


fli) = (4.11) 
and hence 
S(k) = & C, 


where C;, is given by (3.11). The inversion formula 
for the finite transform introduced above is 


Sli) = C, v,{7) (4.12) 


Moreover, the finite transform has the property 
T (pli) f (i)] = — (i) ali) (4.18) 


This could have been used to solve (3.4). 

To return to the continuous case, suppose that 
%(p) is defined in the interval 0 < p < a, where 
a is arbitrary, and, in general, large, and that 
t(p)=0, p>a. Let the interval (0, a) be 
divided into m equal parts, then (4.6) may be 
approximated by 


X,,. = Rpli) X,( pli) 
+ Pi) X,,(p(j)) Ap (4.14) 


where Ap = a/m and p(i) = iAp. This however 
is but Eq. (3.4). It is desired to obtain the solution 
of (4.14) as m +o and it is hoped that this 
will represent the solution of (4.6). To begin, a 
new function p(i) = p(p(i)) is defined 


pli) Ap = 2, — pli) (4.15) 


where A, is that root of the characteristic equation 
(3.7) which satisfies the relation, p(i) << A; < p(i-1). 
Therefore 

0< pi)< 1, 


and the characteristic equation becomes 
5 _pR 4.16 
jar p(k) +k 


remembering that the mole fraction of the jth 
component is z9(p(j)) Ap. If now y is a small but 
otherwise arbitrary number and if « Ap = p, then 
(4.16) can be transformed into 


Pj) 
R= 
Me 


ol J) 


where p(k) zo(k) is some intermediate value of 
p(j) in the intervalk —e Sj Sk+e. As 
m + o, for any A, 0 < A< a, this reduces to 


Anpreas Acrivos and Neat R. AMuNDsSON 


Aag(A) cot pla) 
since 
lim im = cot p(A) 


If finally « is made arbitrarily small and a is very 
large, then formally 


R= p dd + AaolA) w cot wp(A) (417) 


provided that m is always chosen such that 
pm > 1, which is possible. The P associated with 
the improper integral denotes the Cauchy principal 
value of that integral. Equation (4.17) which can 
also be obtained rigorously by using residue theory 
then defines a new function p(A), 0 << A < o, and 
such that 0 < p(A) < 1, which, in a sense, takes 
the place in the continuous problem of the eigen- 
values in the discrete problem. We perceive from 
(4.15) that as m-—> o the discrete spectrum of 
eigenvalues approaches a continuous spectrum 
and as a consequence it would be expected that 
the finite transform (4.10) would pass to an 
integral transform. 

If it is supposed then that z,(p) is a continuous, 
single valued, and integrable function and that it 
does not vanish on the p-axis, then the integral 
transform of a similarly restricted function f(p) 
may be defined as 


T [f(p)] =J0) =p dp 
(4.18) 
where 
w(A) = R — pf ap (4.19) 
One can readily show ne 
T [pf (p)} = Af) — dp (4.20) 


which is the limiting form of (4.18). The inversion 


formula may be formally derived by repeating the 
steps to obtain the characteristic equation, and 


thus 
3 70a) Cd) 
= =pPi ———d 
f(p) = T* [Ff] 
S(p) C(p) 
w(p) (4.21) 


To obtain an expression for the function C(A) 
consider once more the finite problem. For the 
finite problem C, was defined by (3.11) 


Zoli) 
— pli)? 
On the other hand a comparison between (4.21) 
and (4.12) reveals that as the finite problem 
passes to the continuous problem, i.e., AA ~ 0, 


then C, ~ C(A) AA where kA A= A. However, 
as AA is made arbitrarily small 
where p(A) is given by (4.15) and (4.17). But 
= esc? p(A 4.23 
and therefore since arc tan xr = are sin vi = 2 
= (4.24) 


for all positive A. Equation (4.24) may also be 
obtained by more refined methods [10]. 

Equation (4.6) may now be solved for, from 
(4.18), 


T [X,(p)) = X,() = P f x oP) dp 


(4.25) 
it follows from (4.6), (4.20) and (4.18) that 
= RA X,,(A) 
X,(A) = 1 


Therefore 
X,(a) = (RA 
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and from (4.21) and (4.24) 


1 
X,(P) 
(Rp)? 

(Pp) 
which is the solution of (4.6) and where C(p) is 
defined by (4.24). It may be shown that (4.26) 
satisfies (4.14). 

Consider now the case where the function 2,(p) 
is continuous, single valued, and integrable, but 
vanishes only and everywhere in a set of intervals 
BL < p< a,, k =1,2,--¢, where a, > Aye, 
and b, > 6, ,,, and let the remaining part of the 
p axis be the region A. Then the characteristic 
equation is given by (4.17) for A in A but in 
addition one can find ¢ number A, A, in B, such 
that 


+ (4.26) 


=R, 1sk (4.27 
On — p) » lskse (4.27) 
where this integral is proper. Therefore, if 
dp =1 4.28 
Pp (4.28) 


then the solution to (4.6) is 


[a 


(Rp)*"* C(p) (RAY 09 
where X,(p) is defined only for p in A. Finally 
the limiting form of (3.13) is 

a,(p) = . 


1 
(4.80) 
for p in A. 


2. Semi-Continuous Mixtures. Mixtures will 


be called semi-continuous if in addition to con- 
taining an infinite number of components they 
contain a finite number m of components in finite 
concentrations. This is the most general mixture 
that can be considered. Let 2,(i) be the mole 
fraction of the ith distinguishable component 
which is present in finite amount. Then, following 
Bowman [6], the distribution function of the 
overhead product z,(p) may be defined as 


= adp) +2 — p,) di) 


where 5(p — p,) is the Dirac delta function and 
Z4(p) is the distribution function of the continuous 
part of the mixture. Of course 


| ado dp = io dp + = zi) = 1. 


The case will be considered in which 2,(p) is 
continuous, single valued, and integrable along 
the whole p axis and it vanishes everywhere in 
a set of c intervals B as previously defined and is 
non-zero elsewhere, i.e. on the set A. Suppose 
also that p,,1 < m,,isin A and p;,1 <i <m, 
is in B, where m, + m, =m. Then 


= 2p) + — 
+ — pi )ayli) (481) 


and the results of the previous section can be 
applied to the general problem. For, if 


w(A) = R — dp for Ain A (4:32) 
0 


A ofA) 
CO) = + 


ay - R <k S ¢ + m, and 
A,in B (4.34) 


forAin A (4.33) 


forlsksc+m,andAinB (4.85) 
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Then 
C(A) dA + A,*' C, 
(4.36) 
for p in A and 
J A—p a+ 
z,(p) =P oF 4 (4.86a) 


“EAM 

k=1 

for p = pj, 1 S iS my. Moreover from (4.32), 
(4.33), and (4.36) forp =p, 1 Sis m, 

Cy 

— P 


crm 


k=l 


s 
pep) 
R C(A) da + 
A 


z,(p) = (4.36) 
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Similar results may be obtained for the stripping 
column and useful asymptotic expressions, ana- 
logous to (8.26), may be derived. 


NOTATION 


C, = function defined by Eq. (3.11) 
C(A) = function defined by Eq. (4.24) 
D = moles of overhead product per mole of feed 
i, j = component designating variable discrete case 
L = moles of liquid from nth plate 
n = plate designating variable 
p(i) = reciprocal relative volatility for ith compo- 
nent in rectification section 
R = reflux ratio in rectification section 
v,(i) = kth eigenfunction — defined by Eq. (3.8) 
w(p) = function defined by Eq. (4.19) 
W = moles of residue from column bottom 
z,(i) = mole fraction of ith component in liquid on 
nth plate 
2,(p) = distribution function of liquid on nth plate 
Zoli), Zg(p) = composition of total condensate and reflux 
X,, = function defined by Eqs. (3.2) and (4.6) 


y,(i) = mole fraction of ith component in vapour 
from nth plate 
y,(p) = distribution function of vapour from nth 
plate 
a(i) = relative volatility 


Ay = root of Eq. (3.7) 
pp) = function defined by Eq. (4.15) 
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An analysis of the principles of batch and continuous operation 


A. KLINKENBERG 


N.V. De Bataafsche Petroleum Maatschappij, The Hague. (Royal Dutch/Shell Group) 
Summary—An analysis is presented of the meaning of the words “ batch "’ and “ continuous,” 
mainly with regard to transfer processes (also called diffusional processes) and with an extension 
for chemical reactors. In this analysis the various principles of operating such processes are 
examined. It emerges from this analysis, that there is much diversity in the literature with 
regard to the use of such words as batch, intermittent, continuous, step wise, co-current, cross- 
flow, multiple contact and so on. Proposals are made for a handling of such terms, that is believed 
to be a systematic one. 

Finally an effort is made to tie the various arguments that enter into a comparison between 
batch and continuous processing to various aspects recognised in the above analysis. 


Résumé—L auteur précise la signification qu'il donne aux termes “ batch” et continuous,” 
surtout en ce qui concerne les processus de transfert (appelés aussi processus de diffusion) et par 
+ extension dans le cas de réacteurs chimiques. Dans cette analyse, il examine les divers principes 
d'utilisation de tels processus. I] ressort de cette analyse qu’il y a dans la littérature une grande 
diversité dans l'emploi de mots tels que “ batch,” “ intermittent,” “ continuous,” “ step wise,” 
* co-current,” “* cross-flow,” “* multiple contact,” et ainsi de suite. L'auteur fait des propositions 
pour l'emploi systématique de ces termes. Finalement, il fait un effort pour réunir les divers 
arguments qui entrent dans la comparaison entre “ batch ” et “ continuous " en tenant compte 


des divers aspects mis en évidence dans l’analyse ci-dessus. 


INTRODUCTION 


The chemical engineer is confronted with the 
choice between “batch operation” and “ con- 
tinuous operation” in all or nearly all of his 
physical and chemical processes. 

The ideal should be to lay down hard and fast 
rules by the application of which for any specific 
case an unambiguous conclusion can be obtained. 
It has become apparent that such a goal cannot 
be reached, and that the best one can do is to 
enumerate a number of factors or arguments 
which may arise and should be considered. 

These various arguments have been collected, 
for which reference is made to the 11 articles or 
short statements on continuous processing, pub- 
lished in Chem. Met. Eng. in 1945 [2]. 

Recently a symposium was held in the Nether- 
lands by the Royal Institution of Engineers where, 
after a general introduction by VeRscHoor, some 
speakers gave their views on continuization in 
their industries, viz. AssELBERGS (beet sugar), 
GeERRITSEN (salt, caustic and chlorine) and 


Levison (viscose rayon), and where Kramers 


presented a review on holding times. The pro- 
ceedings have just appeared in Chemisch Weekblad 
1954 50 145. 

Some of the arguments raised can be translated 
in terms of money, but not all. The choice is 
therefore not a matter of economy only, but it also 
depends on experience and good judgment. 

It is proposed to discuss here, why the problem 
of opposing “batch” to “continuous” is so 
complex. A very important cause of this com- 
plexity is the fact that the notions “ batch ” and 
“‘ continuous " both are found to be very vague 
and ill-defined. They are used in the literature 
with various meanings attached to them. Liter- 
ally, the words “ batch” and “ continuous ”’ do 
not have opposite meaning, so that it might be 
asked whether they are mutually exclusive and 
together cover all possibilities. 

Another disconcerting feature is the fact that 
there are often many gradations between the most 
discontinuous and the most continuous version of 
a process. For instance, if two materials are con- 
tacted, one may be flowing and the other at 
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rest. Also, pulsating flow is less ‘ continuous ” 
than steady flow. 

A process, which is continuous at first sight, 
will turn out to have a batch aspect if there is a 
preduct which accumulates, such as : 


scale in boilers, 

carbon on catalysts, 

volatile hydrocarbons in circulating volatile 
solvent in mineral oil extraction plants, 

impurities in circulating non-volatile solvents, 
in extraction plants where the products are 


distilled off from the solvent, 
rags at the interface of a continuous settler, 
over-polymerized product in polymerization 
units. 


An indication of the above difficulties is also 
found in the fact that the expressions “ semi- 
continuous and “ semi-batch”’ are sometimes 
used. Thus, Hunter and Nasu [4] describe 
extraction in a counter-current system of mixers 
and settlers as a “ semi-continuous process” as 
opposed to the “continuous” process in an 
extraction tower. 

Waker, Lewis, McApams and GILLILAND 
[17], page 506, use “* semi-batch ” with reference 
to a percolator, where a gas flows through a bed 
of adsorbent. Obviously “ semi-continuous ”’ is 
quite different from “ semi-batch.” 

These arguments are considered sufficient to 
prove that an analysis of the possible meanings 
of the words “ batch ” and “ continuous ” and an 
analysis of possible alternative expressions such 
as “ intermittent’ (see, e.g. WaLker, Lewis, 
McApams and GILLILAND’s intermittent distilla- 
tion [17], page 551-552, and intermittent drying 
[17], page 622) are prerequisites to the discussion 
of the choice between “ batch ” and “ continuous” 
treatment. 

Further to illustrate the necessity of careful 
analysis one may quote one of the authors in the 
1945 Chem. Met. Eng. reports [2], who traced the 
failure of a continuous process to the labourer, 
who was supposed to fill a salt bin, being fast 
asleep. This continuous process failed on a dis- 
continuous aspect. 


The afore-mentioned choice has to be made in 
all branches of engineering. For the present 
treatise some limitation, however, iscon sidered 
desirable. An analysis will therefore be made for 
the important group of unit operations describable 
as “transfer processes” or “ diffusional pro- 
cesses,” the group including distillation, heat 
transfer, extraction, ete. The results will be 
applied to the equally important group of 
chemical reactors. 

Finally, an effort is made to tie some of the 
usual arguments pro and contra “ batch” and 
continuous to the results of these analyses. 

This analysis is founded upon unpublished work 
by W. J. D. van Duck and the author in 1985 
when a nomenclature was developed to be used 
in describing transfer processes. Use thereof was 
made by A. W. J. Mayer in 1987 [4]. 


INTRODUCTION TO A DISCUSSION OF THE 
PrincipLes or Batcu anno Continuous 
OPERATION AS APPLIED IN 
TRANSFER PROCESSES 


The analysis reveals that the word “ continuous ” 
is used in three different ways, sometimes opposed 
to “in stages,’ sometimes with “ batch” as its 
counterpart. 

One form is continuity with respect to place, the 
two others are continuity with respect to time. 

These various possibilities are most fully 
developed in countercurrent processes. These are 
therefore to be examined first and to be followed 
by other processes in the following order : 


Countercurrent processes, 
Single-stage processes, 
Percolator processes, 
Cross-flow processes. 

It should, however, be mentioned that this list 
only covers the simpler possibilities. For instance, 
if percolators are used in series in cyclic rotation 
the action becomes substantially continuous and 
in countercurrent ({17], page 506). Also, bubble 


plate columns show countercurrent between and 
cross-flow on the trays. 
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ANALYSIS OF 
COUNTERCURRENT PROCESSES 

1. Continuity with respect to place 

Upon examination of the compositions* in a 

countercurrent transfer system it is seen, that 

their variation with place is either gradual (in 
packed columns) or stepwise (in plate columns, in 
systems of mixers and settlers). 

It is proposed to denote such processes by the 
names of “ gradual ’’**) and “ stagewise pro- 
cesses, respectively. 

The difference between these groups of pro- 
cesses has a profound influence on the mathema- 
tics of calculating their performance. 

Gradual processes are described by differential 
equations with place as a (continuous) variable. 
Processes in stages are described by finite differ- 
ence equations which, in fact, are simple algebraic 
equations. In the latter, place is a discontinuous 
variable which assumes only whole number values 
(i.e., stage numbers). 

The results, however, obtainable by the two 
processes are very much the same. Any two- 
component separation is obtainable by both. The 
theoretical argument that more-component 
separations cannot be matched because of differing 
rates of transfer for different components seems to 
be of academic interest only. The difference 
gradual versus “ stages is therefore from the 
practical point of view one of mode of operation 
rather than one of principle. 

The literature reveals that : 

(a) there are cases where the word “ batch ” has 
been used to describe a stage-wise or even a 
one-stage operation, and 
gradual countercurrent processes have thus 
far been normally described as continuous 
or differential countercurreni processes. 

Use of “ batch *’—Suerwoop and Picrorp [12], 

page 892, mention “ the single contact in which 


* The temperature in heat transfer will be omitted in 
the following for the sake of simplicity. 

** Unfortunately, the Latin “ gradus " means step. How- 
ever, the development of the meaning of the word “ grad- 
ual,” originally “ step by step ” has been via “ by a great 
number of very small steps ” toward its present meaning 
of a change without steps. 
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the solvent and solution are brought together for 
a single batch extraction.” In the accompanying 
figure this is shown to be carried out in a mixer 
and settler with a continuous feed. The essential 
feature clearly is the single contact, i.e. the single 
stage and the fact, whether the treatment is con- 
tinuous or not, is of secondary importance for 
the discussion. 

MaNLey et al. [7] mention “batch or step 
counterflow extractions, consisting of passing oil 
and solvent in counterflow through a number of 
separate mixers and settlers.” 

STRATFORD et al. [13] mention that “a batch 
method or countercurrent treatment "’ are alter- 
natives to be used in extraction. They make 
“batch ” the equivalent of “ single stage.” 

Waker, Lewis, McApams and GILLILAND [17] 
page 465-466, and Fig. 161, compare “ counter- 
current flow ” to “ batch operation.” 

Use of continuous countercurrent *°—TREYBAL 
[15], page 125, puts “ continuous contact ’’ oppo- 
site to “ stage-wise contact.” 

Maniey et al. [7] mention “ continuous 
counterflow extraction in a single apparatus.” 

SHEeRwoop and Picrorp [12], page 394, oppose 
“continuous countercurrent ”’ to “ multiple con- 
tact countercurrent.” 

The conflicts in the nomenclature in this field 
are clearly brought out by quoting in contrast 
thereto Brown [1], page 301, who says : ‘‘counter- 
current multiple contact may be either batch or 
continuous.’ Here, continuous stands for non- 
intermittent (see page 42). 

Use of “ differential countercurrent *—** Con- 
tinuous countercurrent differential contact,” 
according to Perry [10], page 717 = continuous 
gradual countercurrent. In the next paragraph 
the adjective “ true” is used to distinguish this 
process from multi-stage countercurrent. 


2. Continuity with respect to time 

A convenient starting point is the rather obvious 
statement, to be further refined and modified in 
the following, that “ batch” refers to the feed 
being introduced in separate charges and “ con- 
tinuous " to the fact that the process variables 
(flows, compositions, temperature) are time- 
independent (see [2], pages 101 and 118). 
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Examples are the well-known processes of batch 
and continuous rectification, each of which can be 
carried out gradually (in a packed column) or 
stage-wise (in a plate column). 

For stage-wise processes, not involving a vapour 
or gas a further choice of alternatives becomes 
apparent. Stage-wise countercurrent extraction, 
for instance, may be carried out in two ways : 


1. In separatory funnels and in similar large- 
scale equipment. Here, the same equipment 
serves in regular rotation for all functions 
(mixing, settling, transport). 

2. In mixers and settlers with pumps, where the 
units of equipment do not change their 
function. 


The first of these alternatives falls under the 
description of “ batch ” as given above. If, how- 
ever, one regards batch distillation as a typical 
batch process, it is seen that the extraction in 
separatory funnels may be lacking the character- 
istic features of a typical batch process, such as 
the gradually changing compositions and the 
possibility to collect a great many different 
fractions. In fact, the countercurrent extraction 
in separatory funnels is frequently used to simu- 
late the continuous operation of the other 
alternative, mixers and settlers, the latter being 
more difficult to realize on a small scale. 

If the stages in both are ideal, extraction con- 
ditions are the same, and the extractions are 
continued beyond the initial transient period, the 
products of both extractions are identical. 

Notwithstanding the close similarity to a con- 
tinuous process as described above, the extraction 
in separatory funnels is not continuous in the 
sense that the state at a given point does not vary 
with time. 

Obviously the definitions of batch and con- 
tinuous must be amended to make possible a 
* continuous countercurrent extraction in separa- 
tory funnels.”’ This can be done by using different 
expressions for the continuity of the stage and 
for the continuity of the process as a whole. 

In the following the process using separatory 
funnels will be called “ intermittent” and the 
process using mixers and settlers “non-inter- 
mittent.” The essential feature to distinguish 
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batch and continuous will be the fact whether the 
products of the process are changing or constant. 
In this manner the process in separatory funnels — 
after its transient period — will be a continuous 
countercurrent process in intermittent stages. 

Attention is drawn to the fact that the express- 
ions “ intermittent ” and “ batch ” are both used 
by Waker, Lewis, McApams and GILLiLanp 
[17]. 

In the present treatise usually the names are 
just reversed, as is seen from the following survey : 


W., L., MeA., G. Author 


p. 551 intermittent batch distillation 
distillation 
p. 622 intermittentdrying batch drying 
p. 517 batch distillation cross-flow batch 
distillation 
p. 509 ideal single-stage intermittent single 


batch adsorption stage adsorption 

Our present choice may be justified as follows : 

In the intermittent extraction in separatory 
funnels it is essential that each funnel is used a 
great many times in the same process, i.e. with 
intermissions. In the processes commonly called 
“batch distillation,” “ batch drying,”’' etc., the 
method consists of the working up of a single 
charge. This need not be ever repeated and 
certainly need not be linked up with its inde- 
pendent repetitions, as the expressions “ inter- 
mittent distillation ” and “ intermittent drying ”’ 
would suggest. 

For the simple gradual processes as discussed 
here the distinction intermittent /non-intermittent 
does not arise. There is only batch and continuous 
gradual countercurrent. 


The definitions for continuity do not apply to 
the approach to the steady state of a counter- 
current process. This approach is itself a batch 
process (transient from empty to full as con- 
trasted to transient from full to empty in other 
batch processes). It is of course not the intention 
to let the transient interfere with the definition 
of continuous processes. 


Observations on nomenclature in literature 


1. For transient + continuous intermittent 
countercurrent treatments the name “pseudo- 
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countercurrent ”’ is sometimes used (HUNTER 
and Nasu [5], page 839). 

The “ pseudo * would seem superfluous. A 
certain amount of interpretation to recognize 
countercurrent is always necessary, even in a 
plate fractionation column where descending 
liquid and ascending vapour move through 
different regions of space. 

2. A batch intermittent countercurrent treat- 
ment with a single introduction of feed is 
described by the same authors as multiple 
fractional distribution. 


8. “ Batch countercurrent multi-stage extract- 
tion” (Perry [10], page 718) = continuous 
intermittent countercurrent extraction, i.e. 
“batch ” serves to indicate the intermittent 
operation of the stages. 


4. “Countercurrent multiple contact extraction” 
(= countercurrent extraction in stages) can 
according to SHeErwoop and Picrorp [12], 
page 394, be carried out in batch or con- 
tinuous, exemplified by systems of mixers 
and settlers, as used in batch operation, on 
the one hand, or by the perforated plate 
column on the other. 


These authors add: “At first glance the latter 
would appear to be an example of continuous 
countercurrent flow but in reality should be 
classified as continuous multiple contact, since 
each plate acts as a combined mixer and settler.” 

In our classification the operation of the mixers 
and settlers and the perforated plate column 
would fall under the same heading (i.e. counter- 
current in non-intermittent stages). 


ANALYSIS OF 
SINGLE-STAGE 
For the single stage, consideration of the three 
classification principles as developed for counter- 
current shows : 


1. gradual /stages the process is already lim- 
ited to one stage 
2. non-intermittent/- | identical for 
intermittent 


8. continuous/batch possess ome 
There are therefore two modes of execution of 


a single-stage process (for extraction: mixer + 
settler, and separatory funnel). 

The names non-intermittent and intermittent 
will be given preference because the discontinuous 
process has not the property to make products 
of varying composition. However, continuous 
and batch are in common use. 

Unfortunately, “batch” is often used to 
describe laboratory extractions where the charac- 
teristic feature is “ single stage.’ Authors will 
say that they compare (laboratory) batch ex- 
tractions with (plant) countercurrent extraction. 


Observations on nomenclature in literature 


Single contact is an acceptable equivalent to 
single stage. 

“ Batch treatment” (see SrratrorD et al., as 
cited on page 41) is unnecessarily limited. ‘“* Sin- 
gle-stage batch treatment ”’ is acceptable if it is 
wanted to convey the intention that ‘‘ continuous” 
i.e. non-intermittent operation is not included. 

Simple distillation’ by a “ single vaporiza- 
tion used to separate partially two or more 
volatile liquids ” (WALKER et al. [17], page 516, 
see also [11], page 107) is a term meant to cover 
both single-stage distillation (‘* continuous simple 
distillation,’ such as flashing) and gradual cross- 
stream distillation (“ batch” or “* differential ” 
distillation), the latter to be discussed later. 

This might convey the impression that these 
processes are the continuous and discontinuous 
version of the same, which, of course they are not. 

Discontinuous (i.e. intermittent) single-stage 
distillation is, however, conceivable but it is 
rather impractical because of the large volume 
of a vapour. Its analogue in extraction (the 
separatory funnel) is, however, very common. 


ANALYSIS OF 
PERCOLATION PROCESSES 


In percolators as used in chromatography and 
related processes each phase is contacted suc- 
cessively with portions of the other, as in counter- 
current, but only one of the phases is moving and 
the other is stationary. For this reason the 
percolation in a fixed bed is of necessity a batch 
process. 

THIELE [14] used the name “ percolation ” for 
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processes where a fluid (gas or liquid) is passed 
through a bed of a granular solid. He had in 
mind fixed beds only and there is indeed a need 
for a term describing such operations. WALKER 
et al. ({17], page 465) also refer to percolation 
filters being fixed beds. 

Percolation as defined by Turee is a gradual 
process, but the principle may be extended to 
cover stage-wise processes and these may be 
carried out in non-intermittent and intermittent 
manner. Thus Craic’s well-known extraction 
process [20] is stage-wise and intermittent, and 
the passage of a gaseous mixture through a series 
of washing bottles filled with an absorbent liquid 
is a stage-wise and non-intermittent percolation 
process. 


ANALYSIS OF CROSS-FLOW 
(-STREAM OR -CURRENT) PROCESSES 


The cross-flow principle is easiest recognized, when 
the phase to be treated, say by extraction, is 
flowing continuously through a number of mixers 
and settlers, in each pair of which it is meeting a 
fresh stream of the other phase. This cross-flow 
process is described as stage-wise, non-intermit- 
tent and continuous. By replacing each mixer +- 
settler by a separatory funnel the stage-wise, 
intermittent continuous version is obtained. There 
is also a gradual cross-stream process, where the 
phase to be treated is, on its way, differentially 
extracted, vaporized or the like. 

Upon comparing countercurrent and percola- 
tion, it is seen that the relative movement of the 
phases was the same, whilst the movement with 
regard to the vessel determined whether the 
process could be a continuous process (both 
moving) or had to be a batch process (one at 
standstill). 

These same possibilities exist here : 

In cross-flow as described above, the phase to 
be treated moved from one mixer-settler com- 
bination or one separatory funnel to the next, 
and the process was necessarily continuous. There 
is, however, no necessity to move the phase to a 
new vessel before continuing the treatment. 

By leaving this phase in the same vessel the 
corresponding batch processes are obtained. How- 
ever, all products being the same in the batch and 


continuous version, it is practical to employ the 
name cross-flow for both. The treatment may 
be simultaneously at different places or successive- 
ly at the same place, i.e. place in one case is 
equivalent to time in the other. Gradual cross-flow 
may therefore be continuous in place or in time. 

Thus an ASTM distillation —- when neglecting 
some slight rectification effect — will be termed a 
gradual cross-flow process and the so-called 
multiple-batch extraction is cross-stream in 
stages. For such gradual cross-flow to occur it 
is not necessarily required that one phase is 
removed from the treating vessel. 

In those batch crystallization processes, where 
mixed crystals are formed, the once formed 
crystals are withdrawn from further reaction 
with the liquid by newly deposited layers of 
solid. The continuous crystallizer on the other 
hand operates in one stage and produces a homo- 
geneous product. 


Observations on nomenclature of cross-flow processes 

in literature 

1. The name “co-current” for cross-flow is 
considered ill-chosen ({15], page 129 and 
146; [4], page 340; [6], page 352; [16].) 

2. “Simple multi-stage contact” ([10], page 
716), and “simple multiple contact” ({1], 
page 301) are cross-flow in stages. 


3. “Multiple batch extraction” is also cross- 
flow in stages. 


4. The gradual, cross-flow batch distillation is 
often referred to as differential distillation 
(see [17], page 582), its extraction analogue 
as differential extraction (see [15], page 129). 


5. Waker et al. ({17], page 517), use the ex- 
pression differential distillation, as discussed 
on page 43. 

6. Craic and Craic [20], page 246, describe 
laboratory extractors for the extraction of a 
batch of liquid with a solvent which is 
continuously being recovered by distillation 
(as in the Soxhlet apparatus). Such extrac- 
tion would be a gradual cross-flow batch 
extraction in our nomenclature. It is classed 
as “continuous single stage.” Apparently 
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“ continuous stands for ‘ gradual,’ the “ sin- 
gle stage ” refers to the single vessel, even if 
this is used for a gradual process. 


7. SHeRwoop and Picrorp describe ; “ Multi- 
ple-contact extraction, using fresh solvent in 
each contact "’ ([12], page 392-3). The name 
is clear (= our cross-flow in stages, operated 
continuously). 

8. The use of “ multiple-contact,” without any 
further restriction to describe stagewise 
cross-flow (see [5], page 837) is open to 
criticism. Multiple-contact in fact also 
occurs in stagewise countercurrent equip- 
ment. Thus, TreyBaw (see [15], pages 129, 
146 and 156) rightly distinguishes “ co-current 
(our crossflow) multiple contact’’ and 
“countercurrent multiple contact.” His 
statement (page 129) that countercurrent 
multiple contact is necessarily ‘* continuous ” 
but may be simulated batch obviously 
refers to non-intermittent and intermittent 
treatment, whereas the extraction analogue 
of batch distillation in a plate column has not 
been considered.* 


The following table gives a survey of the 
classification principles discussed so far. 


Subdivision of transfer processes 


-intermittent —-— 
-stages | batch 
Countercurrent -non-intermittent — |— 
uous 
intermittent 
Percolator- non-intermittent always in batch 


(= batch) 
non-intermittent (= continuous) 


batch 
Cin place = continuous 


* Note added in proof : 
Batch extraction in plate colums has been discussed by 
Pratt, H. R. C. ; Chem. Eng. Sci. 1954 3 189. 
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Further criteria of subdivision not to be dis- 
cussed here deal with interchanging the phases 
and choosing the place of the feed, so that 
principles such as rectifying, stripping and com- 
plete rectification, and the analogous concepts 
for cross-flow can be introduced. 


APPLICATION TO CHEMICAL REACTIONS 


A first effort has been made to have some of the 
ideas developed for diffusional processes applied 
to chemical reactions. There is a number of 
chemical reactions between two phases which 
proceed towards an equilibrium just as most 
diffusional processes, and in these cases the 
systematic subdivision can be the same as for 
diffusional operations. 

Here countercurrent and percolation are en- 
countered as important modes of operation. 
Examples are : 


absorption of H,S and CO, in weakly alkaline 
solvents such as K,PO, solution ; 


extraction of weak acids such as mercaptans 
from gasolines with ordinary or solutized 
caustic solutions ; 


simultaneous esterification and distillation ; 
various ion exchange processes. 


By analogy with diffusional processes the above- 
mentioned countercurrent or percolation need not 
be employed if the degree of reaction in a single 
stage is good enough. 

Another important group is found in the 
homogeneous reactions. These may be carried 
out in a single stage (one stirred tank), in more 
stages (series of stirred tanks) or gradually (tur- 
bulent flow in tube). 

In the single-stage process, by analogy with 
transfer processes, the concepts of intermittent /- 
non-intermittent and batch/continuous should 
merge. With the transfer processes the names 
intermittent/non-intermittent were preferred, 
since the two forms can yield the same products 
and do so when the stage is ideal. The occasion 
to compare equipment sized for intermittent and 
non-intermittent operation did not arise in 
transfer processes, because the necesssity to 
contact, exchange and separate causes the 
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intermittent equipment to be different in type 
from the non-intermittent. 

In homogeneous chemical reactors, the same 
type of stirred reactor may serve for both. Here, 
the reaction rate, which may be very slow, is 
solely determining the size. In that case the 
choice of batch versus continuous operation 
becomes of prime importance for the reactor 
size, or inversely, for given size, for the fraction 
reacted. 

A considerable amount of work has been done 
or. the comparison of continuous and batch 
reactors, see, for example, Piret and co-workers 
[3, 8] and Weser [19] covering single, successive 
and simultaneous reactions, both reversible and 
irreversible, of various orders. 

In the single reaction one is interested in the 
fraction converted only. Here, in the continuous 
reactor all reaction takes place at the final con- 
centrations of the reactants, hence slower than 
in the batch reactor (except for zero-order 
reactions). Consequently the continuous reactor 
has the lower conversion for equal volume or the 
larger volume for equal conversion. The difference 
is most strongly felt for higher-order reactions and 
for high percentages of conversion. This dis- 
advantage of continuous reactors can be offset 
to a certain extent by using more than one 
reactor in series, when the earlier vessels in the 
series have the advantage of a higher reactant 
concentration, 

If there are consecutive and simultaneous 
reactions the yields of desired product calculated 
on converted feed material may also differ. 
MacMvu tn [2] has examined this effect for chlorin- 
ation of benzene. He found that the yield of mono- 
chloro benzene is higher in batch chlorination. 

Another argument is met in copolymerization. 

Due to the fact that the different monomers 
have in general different reaction rates, in batch 
copolymerization polymers are formed of changing 
composition. In continuous copolymerization the 
polymer, being formed in a medium of constant 
composition, is uniform as regards chemical 
nature. These predictions have been verified by 
Watt, and Fiori [8]. 

This state of affairs is analogous to the one 
discussed for crystallization (page 44). 


TECHNICAL ARGUMENTS PRO AND CONTRA 

Batcu AND ConTINUOUS 
In discussing batch versus continuous obviously 
all cases where “ batch” has been used in the 
sense of stage and continuous in the sense 
of “gradual” should be disregarded, however 
interesting the comparison of “ stage ’’ versus 
gradual in itself may be. 

Our comparison is to be related to continuity 
in the time. For stage-wise diffusional processes 
in addition to batch /continuous, as related to the 
process as a whole, the principles of intermittent, - 
non-intermittent, as related to the single stage, 
were distinguished. For chemical reactors such 
a further distinction need not usually be made. 

In continuous non-intermittent and in contin- 
uous gradual processes each place in the equipment 
has its well defined function. However, there must 
be equipment for each aspect of the whole pro- 
cess. Thus, a continuous distillation process to 
separate n components requires n-1 columns. 

In a continuous intermittent process several 
dissimilar functions are combined in each con- 
tactor (e.g. mixing and settling). 

In a batch process similar functions under 
different conditions are also combined in one 
unit (e.g. separating component 1 from the mix- 
ture of components 2....m, separating com- 
ponent 2 from the mixture of components 
8....m, and so on, in a batch distillation). 

Consequently, the batch process needs less 
equipment than the continuous process. 

The batch process therefore being less costly 
to install, it is favoured for small units, where 
equipment costs count relatively heavily, and for 
new processes, where pay-out is uncertain. This 
difference is more important the more functions 
the batch equipment has, e.g. with distillations 
to produce a great many fractions. 

The continuous process, being more complex, 
is more difficult to control manually so that 
control is more likely to be fully automatic, with 
a consequent further increase of capital cost. 

The usual comparison therefore is : 
continuous ; fully automatic control, little labour 
batch : less control equipment, labour re- 

quired for maintaining the sequence 


of operations. 
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If, however, the batch process is to be fully 
automatic as well, its control may well become 
more complicated than the control of the con- 
tinuous unit. 

The batch process on account of its flexibility 
is favoured for new and unknown processes and 
for feeds variable in nature, composition and 
amounts. 

These arguments count very heavily if the 
process serves an analytical purpose only (e.g. 
analytical rectification, chromatography). 

In installing a “continuous” process the 
hidden batch aspects due to accumulating re- 
action products, as discussed earlier must be 
considered, Although some types of accumu- 
lation may reach a steady value, the efliciency 
of the process may have been impaired by that 
time to such an extent, that the arguments even 
then count in choosing between batch and con- 
tinuous processing. 

Batches of gas are usually avoided because of 
the low density of gases. , 

The explosion motor is a notable exception and 
there are some rare cases where chemicals might 


be produced by executing a fast reaction many 
times a second in the same container. 

Heat economy by heat exchange is most easily 
obtained in continuous processes, where all flows 
are continuous and the temperatures constant. 
The use of multiple effect certainly requires 
continuity. 

In discussing such chemical reactors, where the 
rate of reaction was a determining factor for the 
size, the continuous reactor had to be larger than 
the batch reactor. The real difference. is smaller 
than the theoretical one, because the batch 
equipment needs off-time for filling and emptying. 

In chemical reactions where the product is an 
intermediate in a sequence of consecutive re- 
actions the choice between batch and continuous 
will influence the ratio of yields of desired and 
by-products. 

It is clear from the above that the choice 
between batch and continuous operation may well 
be such, that the first small unit for a new process 
is advantageously installed in batch and the 
subsequently erected larger unit is a continuous 
one. 
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Book review 


Ullmanns Enzyklopadie der technischen Chemie, 3rd 
edition, by Dr. Forrsr. 


Vol. 3. Abietinsdure bis Arsenverbindungen, XI + 
851 pp, figs., tabs., 19 « 27¢em, bound. DM 108 —, 


Vol. 4. Arzneimittel bis Calcium, NI 
tabs.. 19 27¢m, bound. DM 108 


+ 842 pp, figs., 


Published by Urban & Schwarzenberg, Miinchen — Berlin, 
1953. 


In the third edition of Ullmanns encyclopedia, volumes 3 
and 4 are the first alphabetical ones to be published. Vol- 
ume I,“ Chemischer Apparatebau und Verfahrenstechnik,” 
was published in 1951 as a preliminary to the alphabetical 
series and treated of such matters as unit operations, 
chemical reactors and thermodynamics. Volume 2, due 
this year, will be entitled “ Betrieb und Laboratorium.” 


It was the original intention of the publishers to finish 
two volumes each year. Until now, the rate of appearance 
has been one volume per year. Although it is normal 
practice to publish the parts of an encyclopedia at intervals, 
the disadvantages inherent in this system seem to be 
particularly noticeable in this work. The reasons for this 
ean be found in the rapid growth of chemical industry, 
and in the fact that there is so much detail in these books. 
It is to be hoped that the publishers will overcome this 
difficulty as soon as possible. 


A comparison with the 2nd edition shows in the first 
place a decrease in the number of main headings. There is 
a tendency to group related subjects under one heading 
and this tendency is stronger than it was in the second 
edition. Correspondingly, there is an increased number 
of “see —references. The editor has thus achieved 
improved readability of this text. 


Moreover the subject matter has been conveniently 
arranged without either keeping too strictly to alpha- 


betical order or presenting to the public a collection of 
monographs. In spite of this, much is still to be said in 
favour of the preparation of interim-indexes, especially 
since modern methods of editing and reproduction have 
greatly reduced their cost. 


The type of information given under the various head- 
ings is the same as in the former edition: physical, 
chemical and physiological properties, preparation and 
manufacture, analysis, applications, economics, ete. 
Iixtensive application of flow-sheets and sketches of 
apparatus used in industry greatly contributes towards 
giving a clear presentation of modern manufacturing 
processes. The tremendous growth of chemical industry 
in the past twenty vears clearly shows in such chapters 
as “ Azetylen,’ “Ammoniak, synthetisches,’ and 
Butylalkohole.”” Literature references among which 
those to patents rank first, are numerous throughout both 
volumes. The greater part of all references is given in the 
text rather than at the ends of sections. In the opinion of 
the reviewer, footnotes instead of scattered references 
would improve the readability. 

A comparison of Ullmans work with Kirk and Othmer’s 
Encyclopedia of chemical technology, reveals in the first 
place a difference in the approach to the subject matter. 
The American work is more narrative in character than 
its German counterpart, in which the individual references 
used in compiling the text stand out more clearly. Also, 
in Kirk-Othmer less attention is paid to foreign processes 
than in Ullmann. Volumes 3 and 4 of the latter work 
contain much more technical information (estimated at 
30°.) and less pure chemistry than the corresponding 
American volumes. 

For most cases preference will -urtainly be given to 
Ullmann. To all who are interested in chemical technology 
this work is heartily recommended. 


D. W. van Krevecen 


